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I i ' Abstract: Recent studies, using the AdS/CFT correspondence, of the radiation produced by 

a decaying system or by an accelerated charge in the A/" = 4 supersymmetric Yang-Mills theory, 
^ ■ led to a striking result: the 'supergravity backreaction', which is supposed to describe the energy 

■ density at infinitely strong coupling, yields exactly the same result as at zero coupling, that is, 
l/^ • it shows no trace of quantum broadening. We argue that this is not a real property of the 

. radiation at strong coupling, but an artifact of the backreaction calculation, which is unable to 

, faithfully capture the space-time distribution of the radiation. This becomes obvious in the case 

I of a decaying system ('virtual photon'), for which the backreaction is tantamount to computing 

I a three-point function in the conformal gauge theory, which is independent of the coupling 

^ ■ since protected by symmetries. Whereas this non-renormalization property is specific to the 

^ ■ conformal AA = 4 SYM theory, we argue that the failure of the three-point function to provide 

a local measurement is in fact generic: it holds in any field theory with non-trivial interactions. 
To properly study a localized distribution, one should rather compute a four-point function, 
as standard in deep inelastic scattering. We substantiate these considerations with studies of 
the radiation produced by the decay of a time-like photon at both weak and strong coupling. 
We show that by computing four-point functions, in perturbation theory at weak coupling and, 
respectively, from Witten diagrams at strong coupling, one can follow the quantum evolution 
and thus demonstrate the broadening of the energy distribution. This broadening is slow when 
the coupling is weak but it proceeds as fast as possible in the limit of a strong coupling. 
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1 Introduction 

One topic which has received much attention over the last few years within the context of the 
gauge/string duality is the space-time distribution of the radiation in the strong coupling limit 
of the = 4 super symmetric Yang-Mills (SYM) theory. Originally motivated by studies of 
strongly coupled plasmas in relation with the energy loss by an energetic parton [1-11], this 
problem turned out to be interesting and intriguing for the vacuum case as well, because of a 
surprising result. AdS/CFT calculations of the radiated energy density at infinitely strong cou- 
pling, using the method of the backreaction within the supergravity approximation to the dual 
string theory, led to results which exhibit the same space-time pattern as in the corresponding 
problems at zero coupling: the radiation appears to propagate at the speed of light, without any 
trace of quantum broadening. Originally identified for the case of the synchrotron radiation by a 
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heavy quark [12], this property has subsequently been shown to extend to more general sources 
of radiation [13-21], like an accelerated heavy quark which follows an arbitrary trajectory or the 
decay of a virtual photon. 

The lack of broadening is surprising in that it contradicts our general expectations for 
a quantum theory of interacting fields and, in particular, the experience that we have with 
perturbative studies at weak, but non-zero, coupling. Indeed, in a gauge field theory like = 4 
SYM, one expects the radiation to involve a superposition of quanta with various virtualities, 
including time-like quanta which propagate at subluminal velocity. With increasing time, such 
quanta will separate from each other and also dissociate into other quanta with lower virtualities, 
leading to a spread in the energy distribution along the direction of motion which increases with 
time. At weak coupling, this evolution is well known to lead to parton cascades, in which the 
original virtuality gets evacuated via successive branchings. The associated spreading of the 
parton distribution turns out to be quite slow, because the rate for branching is proportional to 
the strength of the coupling (say, the 't Hooft coupling A = g^Nc in the case of the A/" = 4 SYM 
theory at large iVc). With increasing coupling, the branching becomes more and more effective, 
and the spreading goes faster and faster. In particular in the strong coupling limit A — t- oo one 
expects this spreading to proceed as fast as possible and to occupy the whole region in space 
and time which is allowed by causality and special relativity. 

The following example, to be discussed at length in this work, should illustrate the situation. 
Consider the decay of a 'heavy' photon (an off-shell photon with time-like virtuality) in its rest 
frame. More precisely, the photon is in a localized state represented by a wave-packet centered 
at t = and a; = and which carries a typical 4-momentum = {Q, 0, 0, 0) + 0{l/a), with a 
the width of the wave-packet, assumed to be large: aQ » 1 (see Sect. 2 for details). The photon 
splits into a pair of electrically-charged, massless, partons ('quarks'), which can subsequently 
evolve via 'colour' interactions, that is, by emitting other 'quarks' and 'gluons'. We shall follow 
this evolution to leading order in the electromagnetic coupling, but by letting the strength A of 
the colour interactions to vary from weak to strong. 

When A — )• 0, there is no further evolution, so the final state consists in two on-shell quarks 
propagating back-to-back (by momentum conservation) at the speed of light. The direction of 
propagation of the two quarks is arbitrary, so if one averages over many events one finds an 
energy distribution in the form of a thin spherical shell^ of essentially zero width which radially 
expands at the speed of light: r = t. More precisely, this energy shell has a small width t — rr-^cr, 
which however can be neglected at large times t ^ a. 

If the coupling is non-zero but weak (A ^ 1), the original quarks will be generally off-shell, 
but their virtualities will typically be much smaller than the respective energies. Hence, the 
quarks will propagate with a large boost factor 7^1 before eventually decaying into massless 
quanta. Their radiation will be collimated within an angle ~ I/7 around their direction of 
propagation, leading to a pair of jets in the final state. After averaging over many events, the 
energy distribution has spherical symmetry and a radial spreading t — r which increases with 
time, because of the virtuality distribution of the quanta within the jets. By the uncertainty 

^In QCD, the average distribution has no spherical symmetry because of the bias introduced by the polarization 
vector of the virtual photon. But in A/" = 4 SYM, the anisotropy exactly cancels between the (adjoint) fermion 
and scalar contributions, so the ensuing distribution is isotropic indeed. 
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principle, it takes a time t ~ xQ/fi'^ to emit a quantum with energy fraction x and virtuahty 
/i^. Then, to leading order in perturbation theory, the radial spreading can be estimated as 



where the logarithm has been generated by integrating over the phase-space for the bremsstrahlung 
of a quantum with virtuality k'^ between /i^ and Q^. (For very large times such that A \n{Qt) > 1, 
the higher order corrections become important and will be estimated in Sect. 5.) This argu- 
ment also shows that the typical virtuality /x of the quanta composing the jets is such that 
Aln(Q^/^^) ~ 1, which at weak coupling implies /i <C Q- This confirms that the typical quanta 
are nearly on-shell and thus propagate quite fast: 7 = xQ/ fi ^ 1. 

Consider now the situation at relatively strong coupling, A > 1. Then the virtual photon 
splits into a pair of quarks whose virtualities are comparable to their energies, fj, ~ xQ. These 
quarks are themselves highly virtual and hence they are slowly moving: 7 ~ 1. They will rapidly 
decay into quanta with similar characteristics. We expect this pattern to repeat itself in the 
subsequent steps of the evolution: at each branching, the energy and virtuality of the parent 
parton are quasi-democratically divided among the offspring quanta, which therefore emerge at 
large angles with respect to the direction of propagation of their parent. For sufficiently large 
times t 3> 1/Q, this evolution leads to a parton distribution characterized by a wide dispersion 
in velocities and angles. For the conformal theory A/" = 4 SYM, we expect this distribution to 
be isotropic event-by-event and to show maximal radial broadening, that is, to uniformly cover 
the whole volume at r < t which is allowed by causality. 

Moving to extremely strong coupling A ^ 1, the situation is a priori more complicated, 
since the concept of partons (elementary quanta representing excitations of the quantum fields 
in the Lagrangian which are point-like and nearly on-shell) is probably not useful anymore: the 
matter distribution produced by the decaying photon is made with collective excitations whose 
composition in terms of elementary quanta can be arbitrarily complicated. Yet, since isotropy 
and maximal broadening are already reached for moderate values of the coupling A ~ 0{1), it 
is natural to expect these features to remain valid when A — t- 00. 

These expectations are indeed supported, at least indirectly, by a series of calculations at 
infinitely strong coupling using AdS/CFT. These include studies of the decay of a virtual photon 
using the ultraviolet /infrared (UV/IR) duality [7], calculations of the associated angular corre- 
lations which demonstrate isotropy [22], studies of the jet fragmentation showing the absence 
of point-like partons [23, 24], and also studies of deep inelastic scattering [25-31] leading to a 
similar conclusion: the partons cannot survive in the wavefunction of a hadron, or in a plasma, 
at strong coupling because they efficiently decay towards smaller and smaller values of x. 

Yet, such previous approaches had not address the issue of the radial, or longitudinal, 
distribution of the radiation. For instance, in the study of angular correlations performed in 
Ref. [22], the radial distribution was explicitly integrated over. Also, most of the other studies 
alluded to above were performed in momentum space. The calculation of the backreaction for 
the synchrotron radiation in Ref. [12] is the first attempt in that sense and, as already mentioned, 
it led to the surprising conclusion about the lack of radial broadening. As also mentioned, this 
conclusion applies to other forms of radiation, including our prototype problem — the energy 
produced by the decay of a virtual photon — , for which the backreaction predicts the same 



t — r ~ 
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space-time distribution as at zero coupling: a thin spherical shell expanding at the speed of 
light with a constant width i — r ~ o". This looks puzzling as it suggests that the situation at 
(infinitely) strong coupling could be closer to that at zero coupling, rather than to that at weak 
or intermediate values of the coupling. However, this is not the case, as we now argue. 

A first indication in that sense comes from the following argument, which refers to the 
radiation produced by the decay of a virtual photon. The SUGRA calculation of the backreaction 
amounts to computing a specific three-point correlation function in the underlying field theory, 
which is protected by symmetries and hence it is independent of the value of the coupling. 
Specifically, this correlator reads {Jq Tqo{x) Jq), where Jq is the operator which creates the virtual 
photon (a time-like wave packet of the electromagnetic current operator; see Sect. 2 for details), 
while Tqq{x) is the energy density operator at the 'measurement' point = {t,x). As well 
known, three-point functions in a conformal field theory are fixed by conformal symmetry and 
the (quantum) dimensions of the relevant operators, up to a constant (function of the coupling). 
For the correlator at hand, the operators Jq and Tqq^x) have no anomalous dimensions and 
the overall normalization is fixed by the conservation of the energy. Accordingly, this three- 
point function is independent of the coupling, as anticipated [32]. This property, that we shall 
explicitly check by comparing the respective predictions of the zero-order perturbation theory 
and of the backreaction, 'explains' the lack of broadening shown by the latter, in the sense 
of relating this result to the symmetries of the underlying CFT. But this also demonstrates 
that the three-point function is unable to capture the quantum evolution responsible for the 
radial broadening, since it fails to do so already at weak coupling, where this evolution is well 
understood in perturbation theory. This makes it clear that this three-point function is not a 
good observable for characterizing the space-time distribution of the radiation. 

To summarize, the lack of broadening predicted by the backreaction is not a true feature 
of the radiation at strong coupling, but merely an artifact of computing an observable which 
is not appropriate for that purpose. This observation rises several questions: (i) what are the 
reasons for this failure of the three-point function, (ii) what is the actual physical content of 
a three-point function like {JqToo{x) Jq), and (iii) what are the observables that one should 
study in order to understand the space-time distribution of the radiation. These are clearly very 
general questions and the answers that we shall provide to them are not necessarily new. (Some 
connections with similar problems in QCD will be later pointed out.) But precisely because they 
are so general, these answers are independent of the non-renormalization property of the three- 
point function alluded to above. Most of them apply to any interacting field theory, conformal 
or not, at either weak or strong coupling. 

Specifically, we shall argue that a three-point function like {Jq Too{x) Jq) is truly a forward 
scattering amplitude : the amplitude that the 'target' state created by Jq (i.e. the decaying sys- 
tem) survive intact after interacting with the localized probe operator Tqq{x). In an interacting 
field theory, this amplitude cannot provide information about the internal structure of the target 
at very large times^ t ^ a. Indeed, the quanta composing the target at such late times are very 
soft, as they are the products of many successive branchings, so they cannot provide the high 

^More precisely, we have in mind times which are sufiiciently large to allow for a well developed evolution; at 
strong coupling, the condition f 2> cr is enough in that sense, as we shall later check, whereas at weak coupling 
we require Xln{Qt) > 1. 
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momentum transfer that would be required by a local measurement. (The typical momenta of 
the quanta in the decaying system keep decreasing with t, as we shall see, so they can become 
arbitrarily small for sufficiently large times. By contrast, the typical momenta A'^ transferred 
by the target to the probe are of order l/cr — the maximal value allowed by energy-momentum 
conservation^ — as clear from the fact that the signal has a small width t — r ~ cr.) This argu- 
ment shows that the narrow signal given by the backreaction cannot be a part of the radiation 
in the decaying system at the time t of 'measurement' (the time argument of Tqq{x)). Rather, 
this signal must have been generated at some early time tint ^ before the target had a chance 
to significantly evolve; at that time, the target was composed with only few and relatively hard 
quanta, with momenta A; ~ Q ^ A. But, clearly, such an early emission gives no information 
about the state of the target at the late time t (except at zero coupling) . 

One expects the disparity between tint and t to be maximal at strong coupling, since in 
that case one needs a very small value for tint in order to minimize the effects of the evolution. 
As we shall see in Sect. 3.2 below, this argument is indeed consistent with the calculation of 
the backreaction in AdS/CFT, provided one makes the natural identification between tint and 
the time at which the gravitational wave in AdSs (the 'backreaction') is emitted by the bulk 
excitation representing the decaying system (a SUGRA vector field). 

The above arguments, which explain the failure of the three-point function as a local mea- 
surement, have some other interesting consequences. First, they suggest what should be the 
simplest observable which allows one to study the space-time distribution of the decay: this 
is a four-point function like {Jq Too{xi)Too{x2) Jg) , in which the momentum A transferred to 
the target by the first insertion roo(x2) of the probe operator is then taken away by the second 
insertion Too(2;i). This makes it possible to probe the target with a good resolution (i.e. a 
relatively large momentum transfer A ) without affecting its properties^. Such a measurement 
gives us informations about the state of the target around the space-time point x = {xi + X2)/2, 
with a resolution fixed by the difference xi — X2. 

The previous discussion also tells us under which circumstances a three-point function can 
still act as a measurement: this is possible provided one gives up any radial (or longitudinal) 
resolution, that is, if one integrates over the radial profile of the distribution to get the total 
energy (or, more generally, the energy radiated per unit solid angle), as done e.g. in Ref. [22]. 
Indeed, an operator like the total energy E = J d'^xToo{t,x) involves arbitrarily soft Fourier 
modes, hence it can measure the target without disturbing it. The result of this particular 
measurement is, of course, a priori known: by energy conservation, (Jg E Jq) = Q, with Q the 
energy of the original photon. Less trivial situations occur in the applications of the backreaction 
method to finite-temperature problems. In such cases, one is typically interested in the energy 
deposition in the plasma by a 'hard probe' (a heavy quark, a gluon, or a virtual photon), as 
measured over relatively large space-time scales Ar, At > 1/T, with T the temperature [7, 8, 33- 

^Energy-momentum conservation implies that a forward amplitude like (jj Too (a;) Jq) can receive contributions 
only from the Fourier modes Too (A) whose momenta A*" are smaller than the uncertainty ~ l/cr in the total energy 
and momentum of the decaying system. 

''One should notice the difference between the four-point function that we propose here and the 7i-point 
functions with n > 3 used in Ref. [22]. The probe operators in Ref. [22] are soft, non-local operators, like the 
total energy radiated per unit solid angle along a given direction n : £{n) = lim, ->oo /q°° dtniT°^{t,rn). Such 
operators do not probe the radial distribution of the radiation, but only its angular correlations. 
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40]. Then the method of the backreaction is again rehable, since 1/T is the largest scale for 
quantum broadening in that case. (Indeed, this is the typical value of the broadening by the 
time when the radiation gets thermalized in the plasma.) 

As anticipated, the previous considerations are quite general and in particular they are 
reminiscent of some of the strategies used to study the hadron structure and interactions in 
QCD. Namely, the three-point function and the four-point function above introduced are very 
similar to the electromagnetic form factors and, respectively, the structure functions for deep 
inelastic scattering (DIS), which can both be viewed as measures of the electric charge distri- 
bution in a nucleon, but on very different resolution scales. A 'form factor' is a matrix element 
like {P'\J^{x)\P) where \P) denotes the proton state with 4-momentum and J^(rc) is the 
electromagnetic current operator. For relatively low momentum transfers |A| < where 
A = P' — P and R is the proton (charge) radius, this form factor, which can be studied via 
low-energy electron-proton scattering, provides a good measurement of the proton radius R. 
But if one is interested on the proton structure on much shorter scales r <C -R, as probed by a 
hard scattering which typically breaks the proton, one should rather compute a matrix element 
like n^'"(A) = / d^xe-'^-^{P\J^'{x)J''{0)\P), where A'^ can now be arbitrarily high. This is a 
forward scattering amplitude which via the optical theorem can be related to the total cross- 
section (or 'structure function') for DIS. The experimental measurement of the latter gives us 
the most direct access to parton distributions on short distances. 

Inspired by the above, we shall use here a similar strategy to investigate the space-time 
distribution of the radiation produced by the decay of the virtual photon: we shall compute 
the four-point function describing the DIS between the decaying system and an electromagnetic 
current with space-like virtuality in a boosted frame where the virtual photon propagates at 
nearly the speed of light. (In the context of a decay, this four-point function is also known as the 
'fragmentation function'.) The boost is useful (at least, at weak coupling) to render the parton 
picture of DIS manifest, but our final conclusions at strong coupling can be easily translated 
to the photon rest frame. These results will confirm and substantiate the picture of quantum 
broadening that we previously exposed. 

In the boosted frame, the decaying system looks like a jet — the matter is concentrated 
within a small solid angle Afi ^ 1/7^ around the longitudinal axis (3:3) and within a compar- 
atively small longitudinal interval AX3 <C t behind the light-cone (xa = t) — for any value of 
the coupling. However, at strong coupling this 'jetty' picture is merely the effect of the boost: 
the respective 'jet' is recognized as the boosted version of a distribution which in the photon 
rest frame looks like a uniformly filled sphere with radius r = t. In the boosted frame, this is 
visible in the fact that the longitudinal width Axa of the distribution increases linearly with t : 
Axa ~ (see Fig. 1). This should be compared to the situation at zero coupling, where 

Ax3 ~ (t/7 (the Lorentz-contracted version of a radial width t — r ~ a in the rest frame), and 
also at weak coupling A ^ 1, where Axa increases very slowly with t, as shown in the second 
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line of the equation below (the all-order generalization of Eq. (1.1)) 




for A = 



Axs = {t- X3) 



max 



for < A < 1 



(1.2) 



t 



for A — 00. 
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The above result at strong coupling (the third line in Eq. (1.2)) can be rephrased in a boost- 
invariant way by referring to the typical virtuality /x of the modes in the decaying system: at 
large times, this decreases as /i ~ 1/t . 

Moreover, our analysis of the four-point function will also show that, at strong coupling, the 
matter is uniformly distributed event-by-event within the region of space occupied by the jet, 
meaning that there are no localized substructures, like partons. Indeed, if one tries to scrutinize 
this matter on longitudinal and transverse scales much smaller than its overall respective sizes, 
Ax^ ~ i/27^ and Ax_|_ ~ t/7, then one finds that the fragmentation function is exponentially 
suppressed: it is proportional to exp{— A_|_t/7}, with A_l the transverse momentum transferred 
by the probe current in DIS. By contrast, at weak coupling the fragmentation function is essen- 
tially independent of A_|_, meaning that partons exist and they are point-like. 

So far, we have not been very explicit about the formalism that we shall use and the specific 
calculations that we shall perform. This will be shortly mentioned below, when presenting the 
structure of the paper, and then discussed in more detail in the appropriate sections. As a 
general strategy, we shall perform all our calculations in the framework of the = 4 SYM 
theory, either by using perturbation theory at weak coupling, or the SUGRA approximation to 
the dual string theory at infinitely strong coupling. In particular, we shall use the technique of 
Witten diagrams to evaluate the four-point function describing the fragmentation of the time- 
like photon at strong coupling. A similar calculation has been previously performed in Ref. [41] , 
but only for light -like kinematics (for the 'probe' currents), corresponding to the production of 
on-shell photons. Here, we shall rather focus on the space-like kinematics, which is better suited 
to measure the internal space-time structure of the decaying system. In this paper, we shall 
not address the issue of the stability of the SUGRA approximation against (longitudinal) string 
fluctuations. It has been argued in Ref. [13] that such fluctuations are potentially large and 
unsuppressed in the infinite coupling limit. However, their effects cannot be properly computed 
by lack of a consistent quantization scheme for the string fluctuations in a curved space-time. 
(The heuristic estimates given in Ref. [13] are plagued with severe ultraviolet divergences.) 

Let us also make some comments on the related problem of the radiation by an accelerated 
heavy quark in the fundamental representation of the colour group. There are clearly some 
differences w.r.t. the problem of the decay — notably the fact that the dual object at strong 
coupling is a Nambu-Goto string, instead of a SUGRA field — but we are confident that our 
main conclusions should apply to this problem as well. Indeed, the conclusions concerning 
the quantum evolution at strong coupling, like the maximal broadening, the absence of jets, 
and the absence of partons or other substructures, are universal properties of the radiation at 
strong coupling, which hold independently of the nature of its source. The fact that the radial 
broadening is not visible in the results of the backreaction is again to be attributed to the inability 
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of this method to faithfully capture the space-time distribution of the radiation. To shed more 
light on this point, it is useful to exhibit the CFT correlator which is implicitly computed 
(in the strong coupling limit) via the backreaction. The operator describing the interactions 
between the massive quark and its comparatively soft radiation in the eikonal approximation is 
the Wilson line U{C), with C the trajectory of the quark. Hence, the result of the backreaction 
is proportional to the following correlator in CFT: 

i-(Tr{Z^t(C)foo(x)^(C)}) , (1.3) 

which is recognized as a generalization of the three-point function (jj Too(x) Jq) in which the 
local operator Jq is replaced by the non-local operator U{C). We implicitly assume here a 
large spatial separation between the trajectory C of the quark and the position x of the probe 
operator. (If C is restricted to some bounded region in space with the largest size R, then we 
assume r = \x\ ^ R.) Unlike for (j| Too{x) Jq), we are not aware of general non-renormalization 
properties^ for the correlator (1.3), but this is not essential for our purpose. All that matters 
is that Eq. (1.3) describes an elastic scattering process in which the radiation generated by 
the heavy quark interacts with the probe operator Tqq(x) without being significantly disturbed. 
Then the arguments previously used for {Jq T()o{x) Jq) can be taken over. Namely, the interaction 
with a localized operator is a relatively hard process, which requires a high momentum transfer 
from the target to the probe. The signal carrying such a high momentum can only be emitted 
by quanta which are in the early stages of their evolution, when they are still hard. Such quanta 
have been freshly emitted by the heavy quark and hence they are located in the vicinity of the 
quark trajectory C. Accordingly, the signal carries no information about the structure of the 
radiation at the comparatively remote 'measurement' point x^. This argument is corroborated 
by the backreaction calculation [14] which shows that the emission time tint (identified, once 
again, as the time at which the gravitational wave in AdSs is emitted by the string) coincides 
with the retardation time tr — t — r in the corresponding classical problem — that is, the time 
at which a signal propagating at the speed of light should be emitted by the source in order to 
reach the measurement point r at time t. 

It is finally interesting to mention that correlation functions similar to Eq. (1.3) are com- 
monly used in perturbative QCD to compute the soft radiation produced by energetic partons 
(represented by the Wilson lines), notably in studies of the shape of a jet (see e.g. Ref. [44]). 
However, in such cases the local operator Tqq{x) is replaced with a non-local one, such as the 
total energy radiated per unit solid angle, which is a soft, acceptable, probe. 

The plan of the paper is as follows: In Sect. 2 we shall introduce some general elements 
of the formalism, like the wave-packets describing the virtual photon and the probe operator 
(in both the rest frame of the decay and in a highly boosted frame), and the three-point and 

^This being said, there is empirical evidence that such a property must hold: the results of the backreaction 
in Refs. [12, 14], which include the case of an arbitrary motion for the heavy quark, coincide with the respective 
results at zero coupling up to the replacement A — >■ 4v^ in the overall factor and up to an additional piece 
at (infinitely) strong coupling, which is however a total time derivative and hence averages out for a periodic 
motion. A similar property at the level of the radiated power has been previously observed in Ref. [42]. Such 
non-renormalization properties for the radiation in A/" = 4 SYM, whose precise origin remains to be understood, 
may be viewed as generalizations of similar properties which are known to hold, by conformal symmetry, in 
Euclidean space-time and for simple Wilson loops, like the circular one (see e.g. [43] and references therein). 
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four-point functions that we shall later use to study the decay. We shall explain in more detail 
why a three-point function is not suitable for a local measurement. Also, we shall describe the 
causality constraints on the space-time distribution of the matter produced by the decay. In 
Sect. 3 we shall present the result of the backreaction for the three-point function at infinitely 
strong coupling. With this occasion, we shall correct the original calculation in Ref. [13] by 
adding one term that has been missed there. We shall emphasize the lack of radial broadening 
of the final result and pinpoint the origin of this property in the process of the calculation. 
We shall attempt a physical interpretation for this result in CFT. We shall also perform the 
Fourier transform of the result to a mixed Fourier representation, which is tantamount to using 
a wave-packet for the probe operator. In Sect. 4 we shall present the calculation of the three- 
point function in = 4 SYM at zero coupling (using the mixed Fourier representation, once 
again) and thus obtain exactly the same result as that of the backreaction at infinitely strong 
coupling. Starting with Sect. 5, we shift our attention towards the four-point function that 
describes the DIS of a virtual 7^-current off the decaying system. We first consider the situation 
at weak coupling but late times, where we rely on a leading logarithmic approximation to resum 
perturbative corrections to all orders in Xln{Qt). This will allow us to derive the result for 
longitudinal broadening shown in the second line in Eq. (1.2) and to demonstrate that weakly- 
coupled partons are point-like. Finally, Sect. 6 contains our main new results in this paper, 
namely the calculation of the four-point function at infinitely strong coupling from Witten 
diagrams. For simplicity, that is, in order to avoid a proliferation of diagrams with complicated 
vertices, we shall restrict ourselves to a toy- version of SUGRA — a scalar field theory with 
trilinear couplings. This reproduces the relevant topologies for the Witten diagrams and thus 
correctly captures the physical information which is important for us here: the support of the 
space-time distribution of the radiated matter. We thus find that this matter is uniformly 
distributed over the whole region in space and time which is allowed by causality. 

2 Preliminaries: observables for decaying states 

As announced in the Introduction, our goal is to study the matter distribution created at large 
times by the decay of an unstable excitation of the A/" = 4 SYM theory. For convenience, we 
choose this excitation to be a time-like photon. We follow the standard strategy for introducing 
electromagnetism in AA = 4 SYM, which consists in gauging one of the U(l) subgroups of the 
global SU(4) 7^-symmetry. Then, the electromagnetic vector potentials Acmix) couple to the 
conserved 7?.-current, J^{x), associated with the generator of that particular U(l) subgroup, via 
the action Sint = J d'^x Aemix)Jf_i{x). 

A photon state with given 4-momentum g^, as represented by a plane-wave^ e*'^'^, will be 
on-shell and stable if it has zero virtuality, = 0, but it will be off-shell and unstable when 
its virtuality is time-like, < 0. (The virtuality is defined as q^ = q^q^ = —q^ + q^.) The 
unstable photon will decay into the quanta of A" = 4 SYM which enter the structure of the 
7?.-current (massless fermion and scalar fields in the adjoint representation of the colour group 
SU(A^r;))- These quanta will be time-like too, as they share the virtuality of the original photon, 
so they will themselves decay into other quanta of AA = 4 SYM (including gluons), which will 

^We use a metric convention with the minus sign for the temporal components; e.g. q ■ x = — q^a::" + q ■ x. 
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then spHt again and again, thus progressively evacuating the original virtuality via successive 
branchings. In a conformal field theory like = 4 SYM this branching process will in principle 
go on for ever. If the coupling is weak, the probability for having many splittings is however 
small and the evolution is slow. Then the evolution can be studied in perturbation theory, as we 
shall discuss in Sects. 4 and 5. But at strong coupling, we expect this evolution to proceed as fast 
as permitted by the energy-momentum conservation together with the uncertainty principle. Its 
study can then be addressed within the framework of the AdS/CFT correspondence, and some 
results in that sense will be presented below, in Sects. 3 and 6. 

To be able to follow the space-time evolution of the decaying system, we need to start with 
a perturbation which is localized in space and time. This is conveniently described by a wave- 
packet (WP). Namely, we shall assume that the time-like photon is created by the following 
operator (for more clarity we shall use a hat to denote quantum operators in the CFT) 

Jg = J d^xA^ix)J^{x), (2.1) 

where the 7^-current operator Jfj.{x) is convoluted with a Gaussian WP Aq{x) which encodes 
the information about the 4-momentum, the polarization, and the space-time localization of 
the original perturbation. 

It is instructive to construct this WP in the rest frame of the photon, but then study it in a 
highly boosted frame. This is useful since a boost with a large Lorentz factor 7^1 renders the 
physical interpretation of the quantum evolution more transparent by enhancing the lifetime of 
the virtual excitations (by Lorentz time dilation). In the photon rest frame, the WP is chosen 
as Aq{x) = e^y^(j)Q{x) where e^^^ with A = 0, ±1 are the three polarization states allowed to a 
time-like photon (the polarization index will be omitted in what follows) and 

<t^Q{x) = AAe-^'?* exp|-^^| , j d^x\<t>Q{x)\^ = 1, (2.2) 

is a normalized WP with central 4-momentum = (Q, 0,0,0), which is localized near the 
origin of space-time (t = r = 0) with an uncertainty a. We assume aQ ^ 1, in such a way that 
the Fourier modes k^^ = (k^, k) included in the WP have a typical energy k^ ^ Q and a typical 
virtuality /c^ ~ ^ -Q"^. One has indeed k° = Q + 0{l/a), k' = 0{l/a), i = l,2, 3. 

Consider now the wave-packet in a boosted frame (the 'laboratory' frame) in which the 
photon propagates along the x^ axis nearly at the speed of light. In this frame, the WP has a 
central 4-momentum = {q^,0,0,q^) with v = q^/q^ ~ 1. It is then convenient to introduce 
light -cone components q± = {q^^q^)/ \/2 , in terms of which q^ = (g+ , g_ , 0_l ) and the virtuality 
can be expressed as Q"^ = q^ — q^ = 2q^q^. We shall also need the boost factor, 

^ = y/l\v2 " ^ ^ ^ " ~ V27Q, (2.3) 

which is very large: 7^1. The boosted version of the WP reads Aq{x) = (pqi^) where'' 

^(x) = ATe-'-—'-- exp {-^ - ^ - ^} , (2.4) 



^To simplify writing, we shall not distinguish between lower and upper light-cone components; e.g. A'^ = 
A+ = (A" + A'^)y/2. Also, we use the same notations for the polarization vectors in the rest frame and in the 
laboratory frame, although the longitudinal polarization is of course affected by the boost. 
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with the various widths related to the width a in the rest frame via the following relations, 

a 

cr+ ~ 27 a, a_ ~ — , o-_L = (2.5) 

27 

which express the Lorentz dilation (contraction) of the WP in the x+ (x_) direction. These 
relations imply the inequalities 

cr+g_ > 1, CJ_g+ > 1, cr_LQ > 1, (2.6) 

which in turn guarantee that k^^ ~ for the typical modes included in the WP. The WP (2.4) 
is normalized to unity in the sense of Eq. (2.2) if we choose jTVp = I / {-k'^ a +a ^a]) . 

In order to study the matter distribution produced by the decaying system at late times, 
we shall compute one-point functions like^ 

£g{x) = (Jt foo(x) Jg) , V,{x) = (Jt f++{x) J,) , J,{x) = (Jt J+(x) Jg) , (2.7) 

and two-point functions of the type 

pf)(xi,X2) = (Jtf++(xi)f++(x2) J,), Ilgixi,X2) = (Jt J+(X1)J+(X2) Jg), (2.8) 

where it is understood that all the time arguments xf are much larger than cj+. Recalling the 
definition (2.1) of the operator Jg which creates the state, it should be clear that a 'one-point 
function' like Sg{x) is truly a i/iree-point function in the CFT, and similarly £g'^\xi,X2) and 
Ilg{xi,X2) are truly /our-point functions. 
The integrated quantities 

Eg = J d^x£q{x+,x-.,x±) Pg = J dx^d^x±'Pg{x+,X-,x±) (2.9) 

represent the total energy and the total (light-cone) longitudinal momentum of the state created 
by the operator Jg, and are a priori known: by energy-momentum conservation, they are the 
same as the respective quantities, go and g+, of the original, time-like, photon. In view of 
this, it might be tempting to interpret the integrands in Eq. (2.9), i.e. £g{x) and Vg{x), as 
the corresponding average densities. But this interpretation would be generally incorrect, as we 
now explain. The correlation functions introduced in Eqs. (2.7) and (2.8) are truly (forward) 
scattering amplitudes, which describe the interaction between a 'probe' (operator insertions like 
T++(x) or J4.(xi) J+(2;2)) and a 'target' (the decaying system created by Jg). In the case of the 
three-point functions, this interaction will generally modify the internal structure of the target 
and thus it cannot represent a fine measurement of this structure at the time of scattering. 
The four-point functions, on the other hand, can be used to define a proper measurement, in 
the following sense: if the space-time coordinates xi and X2 of the two operator insertions are 
sufficiently close to each other, then the quantity Ilq{xi,X2) is a measure of the density of TZ- 
charge squared at the central point {xi + X2)/2 as probed with a resolution scale fixed by the 
difference xi — X2 (and similarly for the other four-point functions). 



^The 'average electric charge density' J'q is included here only for illustration: for the problem at hand, where 
the decay is initiated by a electrically neutral photon, we have = in the conformal A/" = 4 SYM theory. 
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The above considerations, to be developed at length in what follows, show that the notion 
of resolution is central to any quantum measurement. This is best appreciated by working in 
momentum space. Then the resolution is controlled by the 4-momentum A'^ transferred by 
the probe to the target, i.e. the momentum carried by the Fourier modes r+_|_(A) of the probe 
operator T^^{x). For a three-point function like Vq{x), energy-momentum conservation requires 
this transferred momentum A to be smaller than the uncertainty ~ 1/cr in the target momentum. 
(For brevity, we use a to collectively denote any of the widths of the target WP, Eq. (2.4). 
More precisely, the conditions on the 4-momentum A'^ of the probe should read as follows: 
A+ < l/cr_ <C q+, A_ < l/o"-|_ ^ q-, and A_l < l/cr.) Yet, in general, it would be wrong 
to conclude that the quantity Vq{x) can be interpreted as the average longitudinal-momentum 
density at x coarse-grained over a distance a. Indeed, even a relatively soft momentum A ~ l/o" 
is still too hard to be absorbed by the target at some large time x+ ^ (7+ and let the state of 
the latter unchanged (within the limits of the uncertainty principle). This is so because, for 
sufficiently large times, the decaying system contains arbitrarily soft quanta. 

This is most easily seen at weak coupling, where one can explicitly follow the evolution of 
the system via successive branchings. One thus finds that the typical longitudinal momenta, 
^4-, of the partons composing the system keep decreasing with time, as expected for a branching 
picture (see Sect. 5 for details). In order to 'see' such partons, a probe should transfer to 
them a longitudinal momentum A_|_ of the order of their own respective momentum A:+. (If 
A_|_ ^ , there is not enough overlap between the probe and the partons to allow for significant 
interactions. If, on the other hand, A+ ^ A;_|_, the probe cannot discriminate the individual 
partons, but only their collective properties averaged over a distance 6x^ ~ 1/A_|..) Clearly, an 
interaction with A+ ~ A;_|_ will strongly affect the struck parton and hence it cannot contribute 
to an elastic scattering unless the momentum transfer A+ is taken back away by a subsequent 
interaction. This can happen in a measurement represented by a four-point function, like 
Tiq{xi,X2), in which case the momentum A transferred to the target by the first insertion J_|_(A) 
of the probe operator is then taken away by the second insertion^ j+(— A). But this cannot be 
the case for three-point functions like those shown in Eq. (2.7). 

We thus conclude that, in order to measure a local quantity, like a density, one can use 
four-point functions, but not also three-point functions. Yet, the latter can be used to measure 
global properties, like the total energy (2.9) : the respective measurement involves no momentum 
transfer, so it cannot affect the decaying system. In general, such a global measurement contains 
no information about the fine spatial distribution of the energy. In some cases, one can recover 
part of this information by exploiting the symmetries of the problem. For instance, the average 
matter distribution produced by the decaying photon has spherical symmetry in the photon rest 
frame. Accordingly, the energy density per unit solid angle is simply obtained as dE/dVl = Eq/An 
with Eq the total energy in Eq. (2.9). But the radial distribution of the energy depends upon 
the detailed dynamics and cannot be inferred in such a simple way. Similarly, the longitudinal 
distribution of the energy in the laboratory frame, i.e. its dependence upon x- = {t — x^)/^/2, 
cannot be deduced without an explicit calculation. In what follows, we shall present such 

^More generally, the 4-momenta Ai and A2 introduced by the two successive insertions, J+(Ai) and J+(A2), 
can be arbitrary but such that their sum Ai + A2 is at most of order l/cr. Via Fourier transform, this sum 
Ai + A2 is conjugated to the central coordinate {x\ + X2)/2 of the measurement process, whereas the difference 
(Ai — A2)/2 is conjugated to the coordinate separation xi — xi and fixes the resolution. 
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x± 




Figure 1. The region in space-time allowed by causality and special relativity for the matter distribution 
produced by the decaying photon in a highly boosted frame. More precisely, the gray band represents 
the boosted version of the half sphere {r < t;x3 > 0} (the region which, in the rest frame, includes the 
quanta whose velocities have a positive third component {v3 > 0). 

calculations for both three-point and four-point functions, at both weak and strong coupling. 

The above discussion shows the importance of simultaneously controlling the localization of 
the probe and its resolution. This can be done by introducing a corresponding wave-packet, i.e. 
by using smeared versions of the probe operators, defined by analogy with Eq. (2.1); e.g., 



The probe wave-packet tpA{y]T) must explore, with the desired resolution, the whole region 
of space where the decaying system can be located at the time of measurement x'^ = r, with 
r ^ A convenient form for the WP is the following Gaussian 



As usual, the central four-momentum A'^ = (A+,A_, A_l) specifies the space-time resolution 
of the probe, whereas the Gaussian controls its localization. The latter is centered at y+ = r, 
with a temporal width which obeys (T+ ^ r (for the time of measurement to be well defined). 
It is furthermore centered at y_ = and y± = 0, with spatial widths (T„ and a± which are large 
enough for the spatial momenta of the typical Fourier components to have only little spread 
around the respective central values: a^A^ ^ 1 and a"_LA_L ^ 1 (compare to Eq. (2.6)). It 
might be tempting to try and enforce the similar condition cj_j_A_ S> 1 on the minus component 
(the light -cone energy), but it turns out that this is not always possible. Indeed, the time 
variable in Eq. (2.11) takes a typical value y+ = r, which is large. In order to avoid the rapid 
oscillations of the complex exponential e~*^-^+ we shall sometimes need to require A_ to be 
small, A„r < 1. Then the condition (T+ A_ ^ 1 cannot be satisfied simultaneously with (T+ <C r. 
But this is not a serious limitation, since we do not need any other temporal resolution scale 
besides the width fT+. To summarize, the WP (2.11) with the constraints alluded to above 
provides a measurement at time r with spatial resolutions 5x± ~ 1/A_|_ and 6x- ~ 1/A+. 




(2.10) 




(2.11) 
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It is finally convenient, before concluding this section, to anticipate the typical resolution 
scales that we shall need in order to probe the structure of the decaying system. This can be fixed 
by comparison with the maximal (transverse and longitudinal) sizes occupied by the system, 
that we shall now estimate. For simplicity, we start in the rest frame of the time-like photon, 
where the matter produced by its decay is restricted to the sphere r < t, simply by causality. 
When boosting this spherical distribution with a large 7 factor, its transverse size remains 
unchanged, that is, Ax± ~ tfip ~ (We used the fact that the time t in the laboratory 

frame is larger by a factor 7 than the time tup in the rest frame.) As for the longitudinal 
extent Axs, this is subjected to Lorentz contraction, yielding Ax^ ~ tRp/l ~ ^/7^- The fastest 
partons propagate at the speed of light, so they will be located on the light-cone X3 = t (or 
x_ =0). Most of the other partons, which are expected to be time-like and thus have velocities 
smaller than one, will be distributed within a region Ax3 ~ i/7^ behind the light -cone. Hence, 
the matter produced by the decay at light-cone time r will be located within a small solid 
angle AO ~ (ArE_|_/T)^ ~ 1/7^ around the X3 axis and within a (relatively) thin longitudinal 
shell Ax_ ~ r/7^ around a;_ = 0. This region is represented as a grey band in Fig. 1. To be 
able to explore its internal structure, we need a probe with sufficiently large spatial momenta 
A_|_ > 7^/r and Aj_ > 7/r. But the opposite case, with A+ <C 7^/t, is also interesting, since 
then the probe measures the matter distribution integrated over the longitudinal (or radial) axis. 
Prom the previous discussion, we expect a three-point function to be a good measurement (say, 
of the energy) when A_|_ <C 7^/t — in which case it correctly provides the energy density per unit 
transverse area (or per unit solid angle in the photon rest frame) — , but not also in the opposite 
case (A+ > ^"^ /t), where the longitudinal resolution is relatively high. These expectations will 
be confirmed by the subsequent calculations, at both strong and weak coupling. 

3 The three— point function at infinitely strong couphng 

In this section we shall briefly review a recent calculation [13] of the three-point function in- 
troduced in Eq. (2.7) in the = 4 SYM theory at (inflnitely) strong coupling, which uses the 
method of the 'backreaction' within the dual supergravity theory. An alternative method, which 
relies on Witten diagrams for supergravity and is perhaps more straightforward to use for the 
calculation of the four-point functions, will be presented in Sect. 6. 

3.1 Backreaction in supergravity 

Within the AdS/CFT correspondence, a time-like photon decaying in the vacuum of the A/" = 4 
SYM theory with infinitely strong 't Hooft coupling (A = g^Nc — )• 00) is dual to a supergravity 
(SUGRA) vector field A^{x^z) which propagates into the bulk of AdSs and whose boundary 
value A^^{x,Q) at the Minkowski boundary [z = 0) is identified with the classical field Aq{x) 
representing the perturbation on the gauge theory side: A'^{x) = (t>q{x) with the wave-packet 
4>q{x) given in Eq. (2.4). 

Within this franework, the first three-point function in (2.7) (the 'energy density' £q{x)) 
can be determined via a backreaction calculation. This refers to the linear response of the metric 
of AdSs to the small perturbation represented by the bulk excitation induced by the boundary 
WP (2.4). In turn, this bulk excitation can be obtained by propagating the boundary field with 
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the help of the relevant bulk-to-boundary propagator (the Green's function for the Maxwell 
equation in AdSs) : 

= I d^2/Z)^^(x-y,z)e,0,(y) = J ^ e'P- D^^-'ip, z) 4>g{p) , (3.1) 

where D^'^{x — y, z) is the Maxwell propagator in AdSs and in the 'radial' gauge = 0. Here, 
z denotes the radial (or 'fifth') dimension in AdSs and we are using the metric (with L the 
curvature radius of AdSs) 

ds^ = GMNdx^^dx^ = ^ [ - dt^ + dr^ + dz^] , (3.2) 

(with M = ;U or z) in terms of which the Minkowski boundary lies at z = 0, as anticipated. 

The SUGRA field (3.1) will be explicitly constructed in Sect. 6.1 below, from which we 
anticipate here the salient features (see also Ref. [13]). Namely, the bulk excitation is a Gaussian 
WP which propagates in AdSs at the 5-dimensional speed of light, with longitudinal velocity 
equal to v and radial velocity Vz = ^/l — v'^ = 1/7- More precisely, at time-^'^ t > cj+, the center 
of the Gaussian is located at 

t , , 

z = —, a;_L = 0, = vt , (3-3) 

7 

with (roughly) time-independent widths fixed by the original Gaussian (2.4) (see Sect. 6.1 for 
details). The physical meaning of the bulk trajectory (3.3) can be understood with the help 
of the UV/IR correspondence [45, 46] : the penetration z of the WP in the bulk is related to 
the virtuality K = ■\/\k'^\ of the typical quanta composing the decaying WP in the boundary 
gauge theory: z ~ 1/K. (For the situation at hand, these quanta are typically time-like: 
k"^ < 0.) Hence, the fact that z is localized near t/^ means that the decaying system at time 
t ^ (T4. involves quanta with a typical virtuality K{t) ~ 7/t and hence a typical longitudinal 
momentum /c+(t) = ^K{t) ~ By the uncertainty principle, such quanta occupy a region 

with transverse area (Aj;_l)^ ~ (V7)^ and longitudinal extent AX3 ~ t/7^ behind the light-cone 
[xz = t). Note that this is the maximal region allowed by causality and special relativity (cf. 
the discussion towards the end of Sect. 2). This qualitative picture for the decaying system at 
strong coupling will be later substantiated, in Sect. 6, by a proper 'measurement' which involves 
the calculation of a four-point function. On the other hand, this picture is not manifest at the 
level of the three-point function Eq{x), to which we now turn. 

The calculation of the backreaction amounts to solving the linearized Einstein equations for 
the (small) change 5Gmn in the metric of AdSs which is generated by the energy-momentum 
tensor t'^'^^ associated with the bulk excitation. Finally, the three-point function (2.7) is inferred 
from the near boundary behaviour of 5Gmn- Mathematically, this is obtained by propagating 
the metric perturbation from the location {x^,z) of its source in the bulk to the measurement 
point xf^ on the boundary (z = 0), with the help of the retarded bulk-to-boundary propagator. 
Strictly speaking, this calculation will yield a retarded three-point function — the retarded 
version of the Wightman function introduced in Eq. (2.7). But this retarded three-point function 
is precisely the physical response function whose space-time localization we would like to study. 

^"Note that x+ ~ y/2t for space-time points located near the center of the WP. 
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For simplicity, we shall replace the bulk WP by a 4-dimensional (5-function with support 
at the central coordinates shown in Eq. (3.3) : t^^^{x,z) oc 5{xz-vt)5^'^\xi_)5{z-th). This 
means that we probe physics on space-time resolution scales which are soft compared to the 
respective widths of the Gaussian WP, which is indeed sufficient for our purposes here. This 
facilitates the calculation of the backreaction, which in general involves an integral over the 
support of the bulk excitation. The result of this calculation reads (see Ref. [13] and also the 
Appendix A to the present paper for details) 

£,{t, x)=^ e r dtt 6{t' -r'- 2{t - vx,)t) 

TT T Jo 

^2qov^ 3^ 2 _ ^2 _ _ ^^3)/)^ (3.4) 

71" 7 Jo 



where r = \x\ and go is the total energy carried by the original WP (2.4) (and therefore also the 
total energy of the evolving partonic system produced by its decay) . Below we shall denote the 
two terms in Eq. (3.4) as £a and £b, respectively, with £q = £a + £b- In the original calculation 
in Ref. [13], the second term £b has actually been missed, so for completeness we shall explicitly 
derive this term in Appendix A. 

Eq. (3.4) can be understood as follows: at time t, the bulk excitation localized at z = t/"f, 
X3 = vt, and x± = emits a gravitational wave SGmn which propagates through AdSs at 
the respective speed of light up to the measurement point x^^ = (t, x) on the boundary. The 6- 
function in the integrand represents the support of the retarded bulk-to-boundary propagator for 
the Einstein equations in AdSs. Its argument follows from causality together with the condition 
of propagation at the 5D speed of light, for both the bulk excitation and the gravitational wave: 

[t - tf = z"^ + {x3 - xzf + x\ = P + - 2xzvi =^ - r"^ = 2{t - vx3)t . (3.5) 

A physical interpretation for this condition back in the original gauge theory will be proposed 
in Sect. 3.2. 

A priori, Eq. (3.4) involves an integral over all the positive values of t, meaning over all 
the values z = t/j of the radial coordinate of the bulk excitation. However, the presence of the 
external derivatives, in the first term and respectively in the second one, introduces an 
important simplification: it implies that the net result for £q comes exclusively from / = 0, that 
is, from the early time when the bulk excitation had been just emitted and was still localized 
near the boundary (z ~ 0). Indeed, after using the (5-function to integrate over t, one finds 

+ 4^ra^*l<'^-W-')l. (3.0) 

where the G-function enforcing r < t (generated via the condition that / > 0) is the expression 
of causality. In the first term, this B-function is multiplied by the factor (t^ — r^) which is linear 
in r^; so the only way to obtain a non-zero result after acting with is that one of the two 
derivatives act on the 0-function and thus generate a (5-function at t = r. A similar discussion 
applies to the second term, which involves an additional factor of (t^ — r^) inside the square 
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brackets and also an additional external derivative. Combining the two terms, one finds 

= -^^3 <5(t2 - r^) . (3.7) 

This describes a spherical shell of zero width which propagates at the ^-dimensional speed of 
light. Retm'ning to the constraint (3.5) on the emission time /, one sees that a signal which at 
time t is located atr = t has been necessarily generated at t = and hence z = 0, as anticipated. 

Now, as it should be clear from the previous discussion, these extremely sharp localization 
properties — the fact that the signal is strictly light-like (r = t) and the (related) fact that 
the whole contribution to the backreaction comes from z = — are to be understood up to 
a smearing on the scale set by the width a of the original WP : in reality, the spherical shell 
has a non-zero radial width t — r ^ a and the values of z contributing to this result are not 
exactly zero, but of order a. Yet, these results — in particular, the fact that the signal appears 
to propagate without broadening (i.e. by preserving a constant radial width up to arbitrarily 
large times) — would be extremely curious if they were to represent the distribution of matter 
produced by a decaying system at strong coupling, as we now explain. 

A thin shell of energy propagating at the speed of light is the result that would be nat- 
urally expected in a non-interacting quantum field theory, or, more generally, to zeroth order 
in perturbation theory for a field theory at weak coupling. In that limit, the time-like photon 
would decay into a pair of (massless) on-shell partons which would then propagate at the speed 
of light. In a given event and in the rest frame of the virtual photon, such a decay yields two 
particles propagating back to back. After averaging over many events, the signal looks like 
a thin spherical shell expanding at the speed of light. In fact, it is straightforward to check 
(and we shall explicitly do that in the next sections) that the result (3.7) of the AdS/CFT 
calculation at infinitely strong coupling is exactly the same as the corresponding prediction of 
the Af = 4: SYM theory at zero coupling. By itself, this 'coincidence' should not be a surprise: 
as explained in the Introduction, the three-point function under consideration cannot receive 
quantum corrections, as it is protected by conformal symmetry and energy conservation. So, the 
corresponding result, as shown in Eq. (3.7), is a priori known to be independent of the coupling. 
But whereas this situation looks natural in view of the underlying conformal symmetry, it might 
still look puzzling from a physical perspective: at non-zero (gauge) coupling, the decay of the 
time-like photon should also involve virtual quanta which propagate slower than light. Then, 
the emerging matter distribution should also have support at points inside the sphere r < t, and 
not only on its (light-like) surface. 

The solution to this puzzle is that, as already argued in Sect. 2 and will be demonstrated via 
explicit calculations in what follows, this three-point function is not a good measurement of the 
energy density produced at late times by the decaying photon. It is clearly a good measurement 
of its total energy, and also of its angular distribution in the photon's rest frame (f = 0), where 
Eq. (3.7) yields the expected result for dE/dil (recall that qq ^ Q in the rest frame) : 



But the radial distribution of the energy density is not correctly represented by Eq. (3.7), in 
any frame. The correct respective distribution will be later computed, at both weak and strong 
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coupling, from four-point functions like those introduced in Eq. (2.8). The results to be thus 
obtained will be very different in the two cases and in particular they will exhibit strong radial 
broadening at strong coupling, in agreement with our general expectations. This being said, it 
would be interesting to understand 'how the conformal symmetry works in practice', meaning 
how is that possible that such a sharply localized result, Eq. (3.7), can emerge from a calculation 
at strong coupling. A possible interpretation for that will be provided in the next subsection. 

For what follows, it will be useful to have a version of the three-point function (3.7) adapted 
to a highly boosted frame (7 S> 1). In that case, it is preferable to work with the probe operator 
T++{x) and the associated three-point function Vq{x), as introduced in Eq. (2.7). At high 
energy, the latter can be estimated as Vq{x) ~ 2£q{x) with £q{x) conveniently rewritten in 
light -cone coordinates. Using 1 — u ~ 1/27^, 1 + f — 2, and hence 

X I (1 — + X_ (1 + f ) X+ r- / ^ 

one finds (with qj^ ~ V^qo) 

2 

n(^) ^ ^ , ""^^ ,3 ^ (2^+^- - 4) • (3.10) 



The denominator in this equation is the reflection of Lorentz contraction, as discussed at the end 
of Sect. 2 : it restricts the longitudinal coordinate x~ to (positive) values satisfying x_ < x+/47^. 
But the presence of the (5-function in Eq. (3.10) entails a much stronger constraint: it implies 
that the signal is localized within an arc of a spherical shell of zero width, or more precisely 
of width (T_ ~ cr/7 (the Lorentz-contr acted version of the respective width in the rest frame). 
This distribution is illustrated in Fig. 2 which should be compared with Fig. 1. One sees that, in 
the boosted frame, the lack of radial broadening mostly manifests itself as a lack of longitudinal 
broadening: the signal (3.10) deviates from the light -cone (a;_ = 0) by a distance x_ = x^/23;_|_ 
(modulo the width o"_ of the shell) which for sufficiently small x± <^ x^/j is much smaller 
than the maximal value ~ x+/7^ permitted by Lorentz contraction. Conversely, this argument 
implies that x± is restricted to values x± < x^/j, which in turn implies that the solid angle 
subtended by the shell is 6il ~ 1/7^- Note finally that Eq. (3.10) yields the correct result for 
the total longitudinal momentum (cf. Eq. (2.9)), as expected: Pq = J dx~d'^x±Vq = q+. 

3.2 A physical interpretation for the 'backreaction' 

As already noticed, the SUGRA results for the three-point function, (3.7) or (3.10), are charac- 
terized by two remarkable and perhaps surprising facts: (i) the signal propagates at the speed 
of light without (radial or longitudinal) broadening, and (ii) the whole contribution to the 
backreaction comes from small values of z < cr. Within the AdS/CFT calculation, these two 
features are related to each other, as we have seen. Namely, the 'backreaction' has support only 
at points satisfying Eq. (3.5), which for small z = t/j < a implies that t — r (or t — X3 in a 
boosted frame) is small as well: t — r < cj in the center-of-mass frame and respectively x^ < (7/7 
in the frame where 7^1. That is, the smallness of z (or of t) implies the propagation of the 
signal at the speed of light. In what follows, we would like to propose a physical interpretation 
for these facts in the CFT. 
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Figure 2. Signal generated from the decaying photon in a highly boosted frame. The width of the grey 
band is (t_ ~ (7/7 (the Lorentz contracted version of the radial width in the rest frame). 

Namely, we shall argue that the interactions responsible for the three-point function are 
highly delocalized in time. The high-momentum transfer between the target and the probe is 
carried by a signal which is emitted by the decaying system at an early time, well before the 
measurement time t at which the signal is absorbed by the probe. This physical emission time 
(denoted as tint in the Introduction) can be identified with the time / at which the gravitational 
wave is emitted by the bulk excitation in the calculation of the backreaction. With this interpre- 
tation, Eq. (3.5) represents the matching condition between the resolution of the probe and the 
kinematics of the target 'partons' which emit the signal. Furthermore, the gravitational wave 
in the 'backreaction' is the AdS dual of the physical signal — a nearly light -like mode with the 
quantum numbers of the probe operator, which propagates at the speed of light from t up to t. 

In order to establish this interpretation, we shall have a new look at Eq. (3.5) which we 
recall is the condition that the gravitational wave propagate at the speed of light in AdSs. For 
a given observation point on the boundary, this condition determines the time t at which the 
gravitational wave is emitted, hence the radial penetration z = t/^ oi the bulk excitation at that 
time and, ultimately, the virtuality K = y^|fc^ of the typical quanta composing the decaying 
system at time t: the UV/IR correspondence implies K ~ \/z = "y/t. For what follows it is 
convenient to fix the transverse coordinate of the observation point — namely, we choose x_l ~ 
(with uncertainty a) — and explore the longitudinal region behind the light-cone (x_ = 0) on 
a resolution scale 5x- which is allowed to vary. As usual, this resolution is controlled by the 
longitudinal momentum of the probe, 5x- ~ and is limited by the longitudinal width 

(T_ ~ (7/7 of the original wave-packet. The best possible resolution 5x- ~ cr_ (corresponding 
to a maximal momentum transfer ~ 7/c") has been implicitly used in the calculation of the 
three-point function 'at a given space-time point', cf. Eq. (3.7) and (3.10). But for the present 
purposes we shall also allow for less precise measurements, with A_|_ <C l/c7_. 

Starting with Eq. (3.5), inserting x_l = and t — = \/2/A_|„, and switching to light- 
cone coordinates, one easily finds (we use the notation r ~ \/2t for the light-cone time of 
measurement and similarly f ~ y/2t for the emission time) 



There are two interesting limiting cases: 

(i) High longitudinal resolution: A_|_ ^ i'^/t. In this case, the probe can discriminate 
longitudinal distances 5x- ~ 1/^+ which are much smaller than the upper limit Ax_ ~ t/j'^ 
enforced by causality and Lorentz contraction. Then Eq. (3.11) implies 

f ~ < r. (3.12) 

^+ 

It is also interesting to estimate (using the UV/IR correspondence) the typical virtuality and 
longitudinal momentum of a quantum from the target at time f : 

Kc.^^-^, k^^V2,K^^. (3.13) 
T 2v27 2 

The last condition (A;+ ~ A_|_) is very interesting: this is the expected matching condition 
between the quanta from the target which emit the relatively hard signal and the resolution 
of the probe. Remarkably, this condition holds here at the 'emission time' f and not at the 
measurement time r. This is in agreement with our expectation that such hard quanta can only 
exist at very early times in the history of the decay. In fact, the above results can be combined 
to yield f ~ kj^/K'^, which is the time interval required via the uncertainty principle for the 
emission of a quantum with longitudinal momentum A;_|_ and virtuality K^. The above discussion 
makes it natural to identify the 'emission time' t (or f) in the SUGRA calculation with the time 
at which the signal measured by the probe operator at time t (or r) has been actually emitted 
by the decaying system in the underlying quantum field theory. 

Returning to Eq. (3.11), let us also consider the other interesting limiting case, namely : 

(ii) Low longitudinal resolution: A_|- <C In this case, the probe cannot discrimi- 
nate any longitudinal substructure, but, interestingly, it can explore the state of the system at 
relatively late times, close to the time r of measurement. Indeed, Eq. (3.11) implies 

2 

f~r, K k+ —. (3.14) 

T r 

This can be understood as follows: the probe is now much softer (A_|_ <C A:+) than the typical 
quanta in the decaying system at time r, so it can interact with the latter without significantly 
disturbing them. 

So far, we have considered a probe with a fixed longitudinal resolution 6x^ ~ that 
is, we have focused on a single Fourier mode, with longitudinal momentum A-|_, of the probe 
operator. But a similar discussion applies to a three-point function in coordinate space, like 
£q{x). The associated Fourier decomposition involves an integral over all values of A_|_, but in 
practice this integral is dominated by its upper limit A_|_ ~ 7/(7, which is the maximal value 
allowed by energy-momentum conservation. Then, Eq. (3.12) implies f ~ and therefore 
z ~ cr. This explains why the whole contribution to the backreaction 'at a given space-time 
point' comes from very small z < a. The fact that the signal propagates at the speed of light 
and without broadening can be qualitatively understood as a consequence of kinematics. Given 
that this signal is carried along by essentially a single mode of the probe — the one with the 
maximal value of A+ — , it naturally preserves a constant width 5x- ~ (7/7. And a signal which 
propagates with t — x^, = const, over a large period of time is necessarily luminal. 
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To summarize, a three-point function with high longitudinal resolution explores the state 
of the target at very early time, much smaller than the time of measurement t. Conversely, the 
only way how a three-point function can measure the state of the target at t is by giving up any 
precision in the longitudinal (or radial) direction. These conclusions will be corroborated by the 
Fourier decomposition of the three-point function to be presented in the next subsection. 

3.3 Momentum space analysis of the backreaction 

In this subsection, we shall compute the Fourier transform of the result in Eq. (3.10) for the 
three-point function in a highly boosted frame. We shall use a mixed Fourier representation 
which involves the component T++(r, A+, A_l) of the probe operator. As explained towards 
the end of Sect. 2, this mixed representation contains the essential information that we need 
about the probe, namely the time of measurement r, assumed to be large (r ^ 0"+), and the 
associated, longitudinal and transverse, resolutions: 5x„ ~ 1/^+ and 5x± ~ 1/A±. 

This change of representation is useful for several purposes. First, it will facilitate the com- 
parison with the zeroth order calculation at weak coupling, to be presented in the next section. 
Second, it will substantiate the argument developed in the previous subsection, concerning the 
correlation between the resolution of the probe and the time of interaction (cf. Eq. (3.11)) . 
Third, it will allow us to explicitly check that the narrow signal seen in coordinate space cor- 
responds to a light-like mode of the probe operator. For the latter purposes, it is preferable to 
perform the Fourier transform before computing the integral over / in Eq. (3.4). 

Consider for illustration the first term, £a, in Eq. (3.4). By using simplifications appropriate 
at high energy, cf. Eq. (3.9), and changing the integration variable from t to z = t/j, we obtain 



^a(t,A+,Ax) = J dx^d^x^e'^+'=—'^^-'^^VA{r,x_,x^) (3.15) 
= ^ J dx.d^x^dzz6"(^2x.{T - V2-fz) - 4 - e^^+"— ^^^-"^ . 

The double prime on the 5-function within the integrand denotes two derivatives w.r.t. its 
argument. It is convenient to rewrite one of them as a derivative w.r.t. z and perform an 
integration by parts to deduce 

£a = 4(z+r I ^ 6' (2x, (r - V2jz) - xi - xM^±^±) e'^+"- • (3.16) 

The Bessel function Jo(A_|_x_l) has been generated by the angular integration over the azimuthal 
angle of x±. We shall now express the remaining derivative of the (5-function as a derivative 
w.r.t. x± and again perform an integration by parts, to obtain (recall that Ji(x) = — dJo/dx) 

= 2,,r f I ' ^ p-^^.MA J,(.))l 

In writing the above, we have also used the (5-function to perform the integral over z in the 
first term within the accolades (the boundary term) and respectively the integral over x± in the 
second term, and we have denoted 



X| fz) = , l2x.T - z(2V2-fX- + . (3.18) 
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The upper limit Zmax in the integral over z is determined by the condition (zn 
yields 

1 T 



0, which 



(3.19) 



Recalling that z = tj^^ ~ t/\/2^ and using < I/A4., this upper limit is clearly consistent 
with our previous estimate for the (maximal) emission time f in Eq. (3.11). 

We now change variables in the integral over z according to z — t- ^ = A_|_X_|_(z), which gives 



dx_e*^+^- 



AxV2x_T 



d^ji(e) 



1 + 



(3.20) 



The Fourier transform of the second term £b in Eq. (3.4) can be similarly computed (in 
particular, this introduces the same upper limit Zmax on z as shown in Eq. (3.19)) and the final 
result reads 



Vq ~ 2[£a+£b) = 8q+ 



7 



/ 



da;_e 



1 + 



(3.21) 



2£g valid at high energy. 



where we have used the relation Vq — 

In order to evaluate the remaining integral over we shall perform approximations ap- 
propriate to the two interesting limiting regimes: A-|- ^ 7^/t and respectively A+ <^ 1^1^- 

(i) Righ longitudinal resolution: A+ ^ 1^ It- In this case, the typical values of X- 
contributing to the integral in Eq. (3.21) obey x_ < 1/A+ ^ t/^"^, so one can neglect the 
second term in the denominator of the integrand. This yields 



' <? — 



7 



dx_e^^+^- Jo(A^v/2^I7) 



tA, 



(3.22) 



The complex exponential can be rewritten as e"*^""^ with A_ = A^/2A+. This relation 
2A-|_A_ = A^ is recognized as the mass-shell condition for a light-like mode. (In fact, if 
one performs the remaining Fourier transform r — )• A_ in Eq. (3.22), one finds a result propor- 
tional to 5(2A_|_A_ — A^).) This light-like mode with high longitudinal resolution is emitted 
at the early time f ~ 7^/A_|_ ^ r and then propagates at the speed of light up to the mea- 
surement time T. The Fourier transform of Eq. (3.22) back to coordinate space is dominated 
by the highest possible values for A_(., namely A™^^ ~ 7/c, which explains why the support 
of the signal in coordinate space lies on the light -cone^^ {x- ~ x^/2x+), with an uncertainty 
~ (T_ introduced by the width of the original wave-packet. 

(ii) Low longitudinal resolution: A+ <C 7^/t". In this case, the typical values of x_ 
contributing to the integral in Eq. (3.21) are determined either by the Bessel function, which 



^^For a generic upper limit A™^", the signal, while propagating at the speed of light, would be shifted from the 
light-cone by a distance 5x- ~ 
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implies < 1/(tA^), or by the denominator of the integrand, which requires < tIi^- 
The last constraint implies that A+3;_ ^ 1 irrespective of the value of A_|_, so we can replace 
giA+x_ ^ The ensuing integral over x_ can be exactly computed by changing variables 
according to x_ = [t lA.^'^)'u} : 



with K2 the modified Bessel function of rank 2. Using [x"^ /2)K2{x) ~ 1 for x ^ 1, we deduce that 
Vq — q+ when A_l <C 7/t. This is as expected: by causality, the decaying sytem has a transverse 
size Ax± ~ r/7 and a longitudinal size Ax_ < ''"/t^ i so when this is probed with much 
poorer, transverse and longitudinal, resolutions, one sees the total energy In the opposite 
limit A_L S> 7/r, the signal is exponentially suppressed (we recall that K2(x) ~ y^7r/2x 
for X ^ 1), meaning that the three-point function does not exhibit any substructure with 
transverse size much smaller than the overall size Ax± ~ r/7. This is again as expected: when 
integrated over x-, the three-point function looks uniform in the transverse plane (at least, at 
points x± <^ t/j) simply by symmetry, that is, as a consequence of the spherical symmetry of 
the signal in the target rest frame. This can be also verified directly in coordinate space: by 
integrating Eq. (3.10) over x_ or, equivalently, by performing the transverse Fourier transform 
in Eq. (3.23), one finds 

n(r,xx) = I dx.V,{x) ^ ^ "'^ (3.24) 

Notice that the low resolution modes are typically space-like : one has indeed A_l ~ 7/r 
and A_ ~ 1/r, hence A^ 2A_|_A_. Consider also the typical values of z and t contributing 
to the signal in Eq. (3.23). By using Eq. (3.19) together with x- ~ t/j'^, one finds Zma.x ~ t/j, 
which implies that / = 72; is commensurable with r. Thus, as already argued in the previous 
subsection, a three-point function with small A+ interacts with the target at times which are 
close to the time of measurement. Yet, because of the low longitudinal resolution, this does 
not bring us any additional information about the state of the system at t. The only physically 
relevant information that we can extract from the three-point function is the energy density 
per unit transverse area, Eq. (3.24), and this is independent of the actual interaction time (as 
it involves an integration over all longitudinal coordinates). 



4 The three— point function at zero couphng 

In this section, we shall calculate the three-point function (2.7) in = 4 SYM in the other 
extreme limit: that of a zero coupling. Our main purpose is to verify that the final result is 
exactly the same as at infinitely strong coupling, as expected from the following facts: (i) in 
a conformal theory like = 4 SYM the general structure of a three-point function is fixed 
by conformal symmetry together with the (quantum) dimensions of the involved operators, and 
(ii) the 7?.-current and the energy-momentum tensor are conserved quantities which are not 
renormalized, that is, they have no anomalous dimensions. Accordingly, the matrix element 
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given in Eq. (2.7) must be independent of the coupling, and this is what shall explicitly check 
in what follows. 

The result of the zeroth order calculation can be easily anticipated. In this limit the time-like 
7^-current decays into a fermion-antifermion (or scalar-antiscalar) pair, which then propagates 
without further evolution. In the center of mass frame of the decay, two back-to-back particles 
moving at the speed of light emerge. The three-point function is not sensitive to correlations 
between the directions of the two decay products, so the answer, in coordinate space, should look 
the same as that of a thin spherical shell of energy whose radius increases with the velocity of 
light. In the boosted frame in which we shall actually do the calculation, the energy distribution 
should be contracted to the part of the spherical shell having solid angle of size 1/7^ around 
the longitudinal axis (the X3 axis along which the decaying current is moving). As we shall see, 
this simple picture is indeed faithfully reflected by the zeroth order result for the three-point 
function. But as we shall later argue, this ability of the three-point function to properly reflect 
the partonic structure of the decaying system is in fact limited to zeroth order: it does not hold 
anymore after including perturbative corrections at weak but non-zero coupling. 

As before, we shall assume that the momentum components. A'*, of the energy-momentum 
tensor T+_|_, are much less than the momentum of the 7^-current initiating the decay. Thus, 
although we are evaluating a transition matrix element, the insertion of affects the decay 
products in such a tiny way that the matrix element corresponds to a faithful determination 
of the average energy flow in the decay. This is of course limited to the present, zeroth order, 
calculation, in which the two partons produced by the original decay do not have the possibility 
to evolve anymore. 

The fact that the three-point function in a conformal field theory is independent of the 
coupling means that, in perturbation theory at least, this quantity cannot correctly describe 
the flow of energy at non-zero coupling, where branchings of the decay products occur. The 
quantum evolution of the partons is on the other hand manifest in the perturbative evaluation 
of the four-point function, to be presented in the next two sections. As we shall see there, this 
evolution leads, at both weak and strong coupling, to the longitudinal broadening of the energy 
flow in the decay. 

4.1 The decay rate 

Our focus in the subsequent calculations at weak coupling will be on the description of the 
average properties of the matter distribution produced by the decay of a time-like 7^-current in 
the N = A SYM theory. To that end it will be useful to have an evaluation of the decay rate V 
of the 7?.-current, an operation which will also allow to introduce our notations. Indeed, in this 
perturbative context, the three-point and four-point functions to be later computed need to be 
divided by T to ensure that they describe properties of a single decay. 

To lowest in perturbation theory, meaning to zeroth order in the gauge coupling g oi J\f = ^ 
SYM and to order in the 'electromagnetic coupling' associated with the 7^-charge, the TZ- 
current can decay into either a fermion-antifermion pair, or into a pair of scalars. To keep the 
presentation as simple as possible, we shall only explicitly evaluate the decay into fermions and 
then simply indicate the changes which occur when adding the scalars. As before, we shall work 
with an 7?.-current boosted along the positive axis, with = Q"^ /2q^ ~ V^^Q and we shall 
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evaluate the rate of decay r+ per unit of light -cone time x^. To the order of interest and for 
the decay into a pair of fermions, this reads 

A,(T,fT 

as illustrated in Fig. 3. The indices A, a, a' refer to the helicities of the decaying 7?.-current, the 
fermion, and the antifermion, respectively. Eq. (4.1) includes a sum over final helicities of the 
fermions and an average (the factor 1/2 in front of the sum symbol) over the initial helicities of 
the current. (The decay rate being the same for any helicity state, we consider here only the two 
transverse helicities: A = ±1.) The phase-space reads d?p = d^p_|_dp+. To evaluate Eq. (4.1) it 
is convenient to use 

d^p dV <5^^^ {q-p-p') = ql z{l - z)dzd(l) , (4.2) 
where = zq^, p'_^ = (1 — z)q^, (p is the azimuthal angle of the fermion, p± = —p'j_, and 

pI = z{1-z)Q\ (4.3) 
One furthermore has 

eu^ip) 7 • e^^^ v,,{p') = e^^h=M= b„„. {a(\ - 2z) - A] . (4.4) 
Using the equations above, one finds 

\ E |en.(p)7-e(")^.'(p')f = 2Q-'[z^ + {I - zf] , (4.5) 

A,o-,o-' 

and therefore 

r+ = /■'dzrz2 + (i_^)2l =^q_. (4.6) 

The decay rate is usually written with respect to the ordinary time variable t in the rest frame 
of the decaying system. Using r_|_rE_|_ = Tt and q-Xj^ ~ one finally obtains 

r = (4.7) 

which is indeed the expected result for the decay of a vector meson with mass Q and purely 
vector coupling of strength e into a pair of massless fermions. 

In A/" = 4 SYM, we also need to include the respective scalar contribution. This is done by 
replacing + (1 — z)^ — 1 in the integrand of Eq. (4.6), so we are finally led to 

rSUSY^f!^^^_ (4.8) 
+ 47r^ 87rg+ ^ ' 

This is the factor which will be used to divide the 3 and four-point functions to get properties 
of the final state normalized to a single decay. 
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Figure 3. Decay of the 7?.-current into a fermion-antifermion pair. 
4.2 The three point function 

We now turn to evaluating the expectation value for the large component of the energy- 
momentum tensor, T_|__|_, at late times in the decay of the time-like 7^-current, to zeroth order 
in the coupling. As in the corresponding calculation at strong coupling, in Sect. 3.1, we shall 
assume that the decay is initiated around the space-time point = 0. In Sect. 3.1, this has 
been enforced by using the wave-packet (2.4). However, as we have seen there, the widths of the 
WP did not play any role in the calculations and in particular they dropped out from the final 
results, like (3.7), because the resolution of the probe was comparatively low (A_|_ < 1/ct_ <^ q+, 
etc). In that respect, the situation will be the same at weak coupling. So, to simplify the discus- 
sion, we shall omit the explicit use of a wave-packet for the incoming 7?--current, but rather use 
its (would-be central) 4-momentum q^^ = {q^,0,0,q^) in order to characterize its localization in 
space and time. 

A similar discussion applies to the probe: strictly speaking, one should use the probe wave- 
packet introduced in Eq. (2.11). But as explained there, the relevant information about the 
resolution and the localization of the probe can be economically taken into account by working 
in the mixed Fourier representation r+_|_(r, A_|_, A_|_). This is precisely the Fourier component 
of the 'backreaction' at strong coupling that we have computed in Sect. 3.3, which will facilitate 
the comparison between the respective results. 

To summarize, in this section we shall compute (with A = (A^,, A_|_)) 

T,(r,A) ^ / d'^e-'^" {Wf+^{r,A)X{0)) 4'^*4'^ (4-9) 

in A/" = 4 SYM at zeroth order in the gauge coupling. The final result of this calculation, after 
being divided by the decay rate F^^^^, Eq. (4.8), will be shown to be identical with the results 
previously obtained at infinitely strong coupling for the quantity ^^(t. A). 

The evaluation of Eq. (4.9) proceeds much as for the decay rate discussed in Sect. 4.1. The 
graph in Fig. 4 shows the energy-momentum tensor interacting with the fermion line, and there 
is a corresponding graph where the momentum A comes off the antifermion line. And there 
are of course also one-loop graphs involving scalar fields to be added at the very end. The 
lines p, p' and p are on-shell, as required by the operator product in Eq. (4.9); this means 
e.g. p- = p\/2p^. This also implies that the 4-momentum A'* exchanged with the probe is 
space-like (A^ > 0), hence the sense of the arrow of time on the corresponding leg is purely 
conventional. (For definiteness, in Fig. 4 we have chosen this line to be outgoing.) 

For the diagram in Fig. 4 (emission from the fermion line), the only differences with respect 
to the calculation given in Sect. 4.1 are a factor of p+ = zg+ (the coupling between the fermion 
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Figure 4. The energy-momentum tensor interacting with the fermion hne. 



and the operator T'+_|_ is proportional to the longitudinal momentum of the former) and the 
replacement of the phase space according to 

d^p d V (27r)^(5(^) {q-p-p') ^ d^p d^p' d^p {2TT)^5^^\q-p-p') 6^^^ (p - p - A) e'^^^" "P" , 

(4.10) 

where within the 3-dimensional 5-function, we have denoted p = (p+,p_l) and similarly for p 
and A. After performing the trivial phase-space integrations using the (5-functions in Eq. (4.10), 
adding the A~emission from the anti-fermion line (this introduces an overall factor of 2) and 
including the corresponding scalar contributions (as before, this amounts to replacing z"^ + (1 — 
z)^ — )• 1 within the integrand), we are left with 

r,(r,A) = dzd^ze-'^P-P~'^\ (4.11) 

where (p is the azimuthal angle between the transverse vectors A_|_ and p± , and 

^ ^ 2p+ 2{p+ - A+) ' ^ ' 

with = ZQ-^ and p± as given in (4.3). Note that A_ = p^ — p-^ is the transfer of light -cone 
energy from the target to the probe. 

So far, we have performed no approximations. At this point we recall that A_|_ <^ q^, 
so unless z is extremely small (which, as we shall see, is generally not the case) we also have 
A_|_ <C p+- Then we can simplify Eq. (4.12) as 

P ± ■ A_L Aj_ p\ 

P+ 2p+ 2p2 



P--P- - — — - 7^A+, (4.13) 



or, after using p± ■ A.± = p±A± cose/), g+ ~ ^/^jQ, and the expression (4.3) for p±, 



p_-p^^\l— (4.14) 



' 1 - z A_L cos A5_ 1 - z a 

2z 7 2z(/_|_ Az 7' 

Inserting this into Eq. (4.11), one can perform the integral over (j) and thus find 



From now on, we shall distinguish, for convenience, between the two kinematical regimes al- 
ready introduced in the discussion of the backreaction: A_|_ j'^/t (low longitudinal resolution) 
and respectively A_|_ » 72/r (high longitudinal resolution). 
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4.2.1 Low longitudinal resolution: A+ <^ 

When A+ <^ 7^/^ and z is not extremely small, both terms in the exponential are much smaller 
than one and hence can be neglected. This is true by assumption for the second term, and it 
is also true for the first term since, as we shall shortly see, at large times one has A_|_ < 7/t. 
(Recall that we consider large times r » 70" > 7/Q.) Then 



which can be exactly integrated (the change of variables z = l/{l + u^) is useful in that respect), 
to finally yield 



As anticipated, the integral over z in Eq. (4.16) is not particularly sensitive to very small values 
2; —7- and the final result in Eq. (4.17) has support at A_l < 7/t. Also, one can check that the 
typical probe kinematics is deeply space-like : the light-cone energy of the probe A_ = p- — p- 
is dominated by the first term in the r.h.s. of Eq. (4.14), which yields A_ ~ A_|_/7 ; hence, for 
A± < 7/r and A+ <C 7^/r, one has indeed A^ > 2A+A_. 

After normalizing by the decay rate (4.8), we obtain the respective quantity for a single 
decay: 



Eq. (4.18) coincides, as expected, with the respective result of the backreaction at infinitely 
strong coupling, presented in Eq. (3.23). As already discussed in that strong-coupling context, 
there is no difficulty in interpreting this result as the average energy measured by a probe 
with strongly space-like kinematics: such a probe has a poor longitudinal resolution, hence it 
measures the energy integrated over the radial profile of the decaying system, within a transverse 
region with radius 5x± ^ 1/A_|_. This energy is correctly given by Eq. (4.18) or (3.23) for any 
value of the coupling. What changes from weak to strong coupling is the radial distribution of 
the energy. In particular, it is only in the zero coupling limit that this radial distribution is 
correctly measured by the three-point function (4.9), as we shall explain in the next subsection. 

4.2.2 High longitudinal resolution: A+ S> 7^/t 

By choosing A_|- ^ 1^ 1^-, one ensures a fine longitudinal resolution in the x_ region populated 
by the decay. To analyze this case, one can again rely on Eq. (4.15), which remains valid so long 
as A_|_ <C q^- Now, however, we cannot neglect the exponential factors in Eq. (4.15) anymore. 
Also, there is no way how the two potentially large phases could compensate with each other, as 
they are both positive definite. So, the only way to avoid strong oscillations is that both phases 
separately remain of order one, or smaller. When applied to the second phase, this condition 
implies that \ — z must be small. By using z ~ 1 together with the change of variables u = 1 — z, 




(4.16) 




(4.17) 




(4.18) 
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one can write 



rr. , .^ e^Q f°° , . f ^^±r\ f Air A+t] 



^ exp z-^-z-^ , (4.19) 



vrA+r [ 2g+ 2A_ 

where the w integration has been extended to n — t- oo because only the small u region is important 
for the integral and we have used the formula 



b 



duJo(aV^)e™'' = le-'^ . (4.20) 



Using A^, ^ it is clear that the dominant phase in the final result in Eq. (4.19) is the second 
phase there. This phase constraints the values of the probe momenta such that (A^/2A+)r ~ 1 
and when this happens the first phase i(A^T/2g+) is much smaller than one and can be ignored. 
For consistency with the previous manipulations, let us notice that the integral in Eq. (4.19) is 
controlled by values of u satisfying 

where the second estimate holds when A_|_/t ~ A^. 

After neglecting the small phase in the second line of Eq. (4.19) and dividing the result by 
the decay rate (4.8), we finally obtain 

r,(r,A) ^ SztV 

rSUSY - A+T ^ > 

Once again, this coincides with the respective result at strong coupling, Eq. (3.22). As explained 
there, the typical value of the light-cone energy (the quantity conjugate to the time of measure- 
ment r) is A_ = A^/2A+, as expected for light-like kinematics. This is indeed consistent with 
our previous estimate A_ = p- — p- for this quantity, as it can be checked by using Eq. (4.14) 
for p_ — p_ together with u = 1 — z from Eq. (4.21). 

What is however specific to the zero-coupling limit at hand is the fact that, in this limit, 
Eq. (4.22) is the Fourier transform of a real measurement. This is possible since the probe is now 
soft as compared to the parton (fermion or scalar) that it interacts with: A+ <^ p+ = zq+ and 
A_|_ ^ p±. (The second condition follows by using Eq. (4.3) and Eq. (4.21) to successively write 
p± ~ ^/uQ ~ (g_|_/A+)A_L » A_L.) So, for this particular problem, even a probe with relatively 
'high resolution' (which can discriminate longitudinal and transverse structures much smaller 
than the maximal respective sizes, Aa;_ ~ and Aa;^ ~ permitted by causality) is 

still soft enough to provide a coarse-grained measurement over a volume much larger than the 
volume occupied by the struck parton. Since 2; ~ 1, the struck parton carries most of the original 
photon energy: p+. = zq+ ~ g+. Accordingly, we expect the result (4.22) to equal this energy 
q+. times the probability for the trajectory of the (small) parton to intersect the (comparatively) 
large area of the 'detector'. And indeed, the prefactor in (4.22) can be given such a simple 
geometric interpretation, as we now argue. 
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Namely, the decay occurs over a solid angle 1/7^. The measurement covers a region with 
transverse area 1/A^ and hence a solid angle ~ l/(A_Lr)^. Thus the measurement covers a 
fraction (7/A_|_r)^ of the solid angle of the decay. When ~ A_|_/r, this is the same as the 
prefactor in (4.22) except for the factor Hence, Eq. (4.22) is the fraction of the energy of 
the decaying system which propagates within the solid angle covered by the 'detector'. 

5 Jet evolution at weak coupling but late time 

With this section we begin the study of the four-point functions introduced in Eq. (2.8), first in 
the context of perturbation theory at weak coupling. As anticipated in Sect. 2, these correlations 
correspond to measurements which can actually probe the space-time distribution of matter 
produced by the decaying system. This will be manifest in the subsequent discussion of the 
situation at weak coupling, where we shall see that the four-point functions reveal the partons 
and their evolution. First, in Sect. 5.1 we shall introduce the partonic fragmentation function 
of a jet and discuss its evolution at weak coupling but late times. Then in Sect. 5.2 we shall 
relate this fragmentation function to a specific four-point function (essentially, the function Hq 
in Eq. (2.8)) that can be measured (at least in principle) via deep inelastic scattering. 

5.1 The general picture 

We now turn to the case where the coupling constant of A/" = 4 SYM is small, but non-zero. 
As before, we are interested in the matter distribution produced by the decay of a time-like TZ- 
current at sufficiently large times — much larger than the characteristic time scale tq = 2q+ / 
for the occurrence of the first decay. If one waits for such long time, the two high-momentum 
partons produced by the original decay must have evolved into a large number of softer partons. 
In general if one wishes to keep track of the number of partons in the evolution, both their 
small-x and their DGLAP evolution are important. However, if one only wishes to follow the 
time-dependence of the energy distribution, the DGLAP evolution is sufficient. Indeed, at weak 
coupling at least, the small-x partons, although more numerous, carry only a negligible fraction 
of the totale energy. 

The DGLAP evolution characterizes the change in the parton distributions (also known 
as 'fragmentation functions' in the context of a time-like evolution) with decreasing virtuality, 
from the original virtuality of the 7?.-current down to the virtuality scale of interest 
This can be also viewed as an evolution in time, by using the relation between the lifetime 
of a parton generation and their virtuality given by the uncertainty principle (see Eq. (5.6) 
below). The fragmentation function D{x,Q'^ / fi'^) represents the number of partons of a given 
species (which for our present purposes can be left unspecified) per unit x and with virtuality 
comprised between and As we shall explain in the next subsection, this quantity truly 
corresponds to a four-point function which can be measured via deep inelastic scattering. 

In what follows, we shall assume that //^ <^ and we shall limit ourselves to the 'leading- 
logarithmic accuracy' (LLA), in which the DGLAP equation resums powers of Aln(Q^//i^) to all 
orders. This equation is most conveniently solved via a Mellin transform x ^ j. This introduces 
the 'anomalous dimension' (the Mellin transform of the DGLAP splitting kernel), which to 
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the accuracy of interest reads [47] 

^0') = (5.1) 

where ■0(7) = d lnr(7) /d7 is the di-gamma function. Then the fragmentation function D(x, Q'^/ jj^) 
is expressed as the inverse Melhn transform 

x^D{x,Q'/fi') = f ^e(^-2)i-{V-)+70-)in{QVM^)^ (5.2) 
J 27r? 

where the j-integration goes parallel to the imaginary axis and to the right of j = 1. We 
used the initial condition D = 6{x — 1) when = fj,^. The factor x'^ has been introduced for 
convenience. Indeed, we are mainly interested in finding what are the x-values of the partons 
which carry most of the energy of the decay. To that aim, it is useful to multiply D{x, Q'^ / fi^) by 
x'^, one factor of x to get an energy weighting and another one to count the number of partons. 
(With the present conventions, the number of partons in an interval dx is dA^ = Ddx.) 

We assume that Q^///^ is very large, such that Aln((5^//U^) ^ 1, and then the integral over 
j can be evaluated in the saddle point approximation. The saddle point js obeys 

lni+7'(j,)ln^ = 0, (5.3) 
and leads to 

x'D{x,Q'/f,') ~ (^^J . (5.4) 

As already mentioned, we are interested in the values of x which maximize (5.4). The maximum 
value of the function /(j) = — [j — 2)^'{j) is /(2) = 0, so the values of x which dominate 
the energy-momentum sum rule is given by Eq. (5.3) with jg = 2, namely 

In - = -7'(2) In ^ = 4 In ^ . (5.5) 

We see that Xc gets smaller as /x^ gets smaller, as expected in view of our physical picture 
of parton branching. Now, since the virtuality decreases along the branching process, it is 
convenient to express /x^ in terms of the duration r of the decay. To that aim, we shall observe 
that, by the uncertainty principle, the partons with longitudinal momentum fraction x and 
virtuality /u^ have a formation time 

Tform - 5- • (5.6) 

To be able to produce the partons with a given x and /U^, the evolution must occur over a time 
T = Tformix, fJ-"^)- (If T ^ Tformi such partons have no time to be formed, whereas if r » Tform) 
then they have already decayed by the time of measurement into partons with smaller values for 
X and /.f^.) These considerations show that, for a given x, one can express the lower end fi'^ of 
the virtuality evolution in terms of the evolution time r, by equating Tfoj.^{x , fi"^) with r. This 
yields the following relation 

- = x% (5.7) 

To /^^ 
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between r, x and /i^. (We recall that tq = 2q^/Q^ is the formation time for the first decay of 
the 7^-current into a pair of partons.) Using this relation for x = Xc, one can finally rewrite 
Eq. (5.5) as an equation for the evolution of Xc with r : 

A/24 

- • (5-8) 

This result explicitly shows which are the partons which carry most the energy of the decaying 
system at a time r after the decay has begun. As expected, Xc decreases with time, albeit only 
slowly (as a small inverse power of r), because the evolution is slow when the coupling is weak. 

5.2 The four point function and deep inelastic scattering 

We now turn to an analysis of the decaying state by performing a deeply inelastic scattering, 
at time r, on that state. We shall use an 7^-current not only to create the decaying system but 
also as a probe to measure this decay via DIS. To better control the space-time resolution and 
localization of the probe, we shall associate a wave-packet to the respective 7^-current, 

Mt,A) = J dSijA{y;r)J^{x), (5.9) 

with the probe wave-packet ipA{y',T) as introduced in Eq. (2.11). To ensure a good resolution, 
this wave-packet needs to be strongly space-like (see below for the precise conditions). The 
'deep inelastic scattering' is the process where the decaying time-like system absorbs the space- 
like probe current and thus evolves into some arbitrary final state. The inclusive cross-section, 
also known as the DIS structure function, is obtained by summing over all the possible final 
states. Via the optical theorem, it can be related to the following forward scattering amplitude, 
which is a Wightman function 

n,(r,A) = / d'xe-'"-^ ( J^(x) J+(r, A) J+(r, -A) 4(0)) 6^*6^. (5.10) 

A 

As manifest from the above writing, the 4-momentum A'^ transferred to the target by the first 
insertion of the probe current is then taken away by the second insertion, so the target can 
be measured with high resolution without being perturbed, as anticipated in Sect. 2. Strictly 
speaking, the above statement refers only to the central value A'^ of the probe WP 4-momentum, 
but in this case one can chose the WP to be strongly peaked in momentum at this central value, 
with negligible spread. 

For reasons to later become clear, it is now preferable to choose A+ = and use the other 
components of the probe momentum, Aj^ and A_, to control the transverse, longitudinal and 
temporal resolutions of the experiment. These non-zero components can be arbitrarily large 
and they have negligible spread, meaning that the corresponding widths are relatively large: 
o"_lA_|_ S> 1 and (7_|_A_ ^ 1. More precisely, we shall chose these widths large enough for the 
detector to cover the whole spatial region where the decaying system can be localized at time 
T, in order not to miss any parton; this requires d'± > t/^ and cj+ > Also, as before, we 

require ^ r in order for the time of measurement to be well defined. As we shall see, there 
is indeed no difficulty to satisfy all these conditions for the problem at hand. In particular, the 
probe is (strongly) space-like, A^ = A^ > 0, as anticipated. 
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Figure 5. The scattering of the current on a parton with momentum k. 

The perturbative analysis of deep inelastic scattering at weak coupling is well developed in 
the literature and will be not repeated here, especially since the corresponding result is already 
known to the accuracy of interest: the structure function (5.10) is proportional to the partonic 
fragmentation function introduced in Sect. 5.1. In what follows, we shall simply explain the 
relation between the kinematics of DIS and the variables x and Q^/^^ of the fragmentation 
function. To that aim, we consider the absorption of the current with momentum A by a 
parton with momentum k, as illustrated in Fig. 5. The parton can be assumed to be on-shell 
both before and after this interaction. This is implicit in our (Wightman) prescription for 
ordering the operators in the four-point function. It is also physically reasonable, since the 
virtuality k"^ ~ fi'^ of the parton is much smaller than its longitudinal momentum squared fc^ 
with fc+ = Hence, we can write k'^ ~ {xq^,0,O±). The on~shellness conditions = and 
{k + A)^ = then imply A^ = 2xq^A-, thus fixing the longitudinal momentum fraction x of 
the struck parton. The corresponding virtuality fi^ can be estimated as in Sect. 5.1, by equating 
the time of measurement r with the formation time (5.6). 

The above considerations motivate the following expression for the four-point function at 
hand: 

^^Ug{T,A) = xD{x,Q^/ii\Al/iJi^), (5.11) 

where 

x= and (5.12) 

2g+A_ T 

The fragmentation function in the r.h.s. of Eq. (5.11) refers to partons with longitudinal fraction 
X, virtuality /i^ and transverse area 1/A^. The normalization factor in the l.h.s. can be 
understood as follows: the fragmentation function counts the number of partons in one decay, 
whereas the four-point function (5.10) provides an integrated version of this quantity over the 
typical duration 6xj^ ~ ^1 ^- of a collision between the probe and the target. Hence, to obtain 
the number of partons per decay, one needs to divide ng(r. A) by the typical number of decays 
occurring during an interval (5x+, namely (5x+r^^^^. 

By the uncertainty principle, a parton with longitudinal momentum = xqj^ is delocalized 
over a distance /S.x_ ~ 1/xqj^. Since our probe can actually 'see' such partons, we deduce that 
it has a longitudinal resolution 5x- ~ 2A_/A^. This can be made arbitrarily small by taking 
Ax to be sufficiently large. In particular, for the typical partons that carry most of the total 
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energy at time r, Eq. (5.8) implies 



1 1 /t\^/24 

Sx^> ^ - - • (5-13) 

Xcq+ q+ \ToJ 

This grows with r, albeit only slowly: as intuitive from the picture of parton branching, there is 
a spreading of the radiation in the longitudinal direction with increasing time, but this spreading 
is slow, since so is the evolution. 

In order for the scattering to give a well defined value for x as indicated above, it is important 
that the temporal uncertainty in the wave-packet, cj+, obey 

r»a+»^:.||t. (5.14) 

Once again, this is easy to achieve so long as A_|_ is not too small. (In particular, one needs 
A_L ^ /i, as clear by comparing Eqs. (5.14) and (5.11).) By the same token, the temporal 
resolution ~ 1/A_ of the scattering is extremely good: the space-like current is absorbed 
over a very short time Sx+ <C o"+, hence it probes the state of the decaying system at r. 

Consider finally the dependence of the fragmentation function upon the parton transverse 
size 6x± ~ 1/A±. This has not been mentioned in Sect 5.1, since it is almost trivial in the 
present context: the weakly-coupled partons are point-like (they occupy a negligible area in the 
transverse space), hence the structure function is independent of the probe transverse resolution 
A_L (for a given value of the longitudinal momentum fraction x). More precisely, there is a weak 
dependence, via powers of Aln(A^//i^), which has been neglected in the discussion in Sect 5.1 
and which is generated by the quantum evolution between the virtuality scales fi^ and A^ of 
the struck parton and of the probe, respectively. (This corresponds to radiative corrections to 
the interaction vertex in Fig. 5 and is controlled by the space-like anomalous dimension of the 
probe.) But such a weak dependence, which reflects the perturbative evolution of the partons, 
does not alter the basic fact that partons are essentially point-like^^ . 

To summarize, perturbation theory at weak coupling but large time predicts that the de- 
caying system evolves via successive parton branchings into partons which are point-like in the 
transverse plane and whose longitudinal spreading (in the sense of their deviation 6x- from the 
light-cone X3 = t, where all the particles would be located in the limit of a zero coupling) is 
slowly increasing with time. 



6 Witten diagrams at strong coupling 

In the previous section we have seen that, at weak coupling, the four-point function (5.10) 
describes the quantum evolution of the decaying system via parton branching, whereas the 
three-point function (4.9) cannot do so. In what follows we would like to extend these findings 
to strong coupling, by showing that the four-point function computed within AdS/CFT from 
Witten diagrams shows indeed quantum broadening and no trace of point-like substructures 

^^If the situation was difTerent, that is, if the partons had some intrinsic, time-dependent, transverse size r(r), 
then the fragmentation function would exhibit a strong, geometric, dependence upon the variable A±r(T) and in 
particular it would rapidly die away when A±r{T) — > 00 (since in that limit, the probe would be unable to see 
the partons anymore). 
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('partons'). To that aim, we shall focus directly on the space-like kinematics for the probe, 
since this is the kinematics which has revealed partons at weak coupling. Besides, we know by 
now that it is only with this kinematics that we can access the state of the system at a time 
close to the measurement time. 

The Witten diagrams involve interactions occurring inside AdSs between the bulk exci- 
tations induced by the 'target' and 'probe' operators from the boundary gauge theory. The 
interaction vertices are local in AdSs and they can be connected via bulk-to-bulk propagators. 
For simplicity, we shall perform our calculations in a scalar model for SUGRA which is a scalar 
theory for a massless field in AdSs with trilinear interactions. That is, the 7^-currents from the 
previous discussion will be replaced (for both the decaying system and the probe) by scalar, 
'dilaton', operators, denoted as O. This model generates Witten diagrams with the same topol- 
ogy as the relevant SUGRA diagrams [41], but of course there are fewer such diagrams and with 
simpler vertices. Although strictly speaking we modify the theory by doing that, we do not 
believe that this could alter our conclusions. Indeed, here we are only interested in very robust, 
qualitative, properties like the support of the four-point function as a function of the probe 
kinematics. For a space-like probe at least, such properties are captured (in the economy of the 
SUGRA calculation) by the various bulk-to-boundary and bulk-to-bulk propagators, but they 
are not sensitive to the detailed structure of the vertices. 

To start with, in Sect. 6.1 we shall present a careful construction, using wave packets, of the 
bulk excitations corresponding to the target and the probe. This will allow us to check some 
approximations used in the previous manipulations, in particular the fact that one can treat 
the bulk excitation associated with the decaying system as a 'particle falling in AdSs'. Then, 
in Sect. 6.2, we shall compute the three- and four-point functions in the scalar model. The 
calculation of the three-point function is shown only for completeness, namely to demonstrate 
that, even for such a scalar toy model, the Witten diagram provides a result which is qualitatively 
consistent with the Fourier transform of the backreaction computed in Sect. 3.3. The calculation 
of the four-point function confirms what we have been so far advocating, that jet evolution at 
strong coupling leads to total, radial and transverse, broadening, with no trace of substructures. 
This will be emphasized in the physical discussion of the results, in Sect. 6.3. 

6.1 Preliminaries: bulk excitations 

As already explained in Sect. 3.1, the bulk excitations representing the decaying system and 
the probe are obtained by propagating the respective boundary fields towards the interior of 
AdSs with the help of boundary-to-bulk propagators. From now on we shall restrict ourselves 
to scalar perturbations, corresponding to 'dilaton' operators in the boundary field theory. 

Consider first the decaying, time-like, system. The corresponding boundary wave-packet 
will be taken as in Eq. (2.4) . The associated bulk excitation reads 



$g(x,z) 
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(6.1) 



(2vr)4 



where (with TV' = TV {2Tr)'^ a+a^a'j_) 




(6.2) 
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is the momentum-space version of Eq. (2.4) and (with = p^p^ = — 2p+p_ 



D{p,z) 



2 ■K2{Pz) if P'>0, 



(6.3) 

H^'^(|P|z) if P^<0, 



4 

is the time-ordered (or Feynman) boundary-to-bulk propagator in momentum space. (For the 
space-like modes and also for the time-like ones with positive energy, this coincides with the 
respective retarded propagator.) Given the conditions (2.6) on the widths of the WP, it is 
clear that the typical momenta allowed in the integral over p^ in Eq. (6.1) are time-like, with 
-p2 ~ g2 ^j^^ ~ > 0. 

As before, we are interested in large times x+ ^ 0"+ ^ Then, as we shall shortly see, 

the bulk excitation is localized at relatively large values of z, such that Qz ^ 1. Accordingly, one 
can use the asymptotic form, valid for \P\z ^ 1, for the Hankel function within the propagator: 



H2^^(^) - A/ — e*'^"'^^''/^^ when x > 1 . (6.4) 



The momentum integral in Eq. (6.1) is then controlled by the condition that several large and 
strongly oscillating phases compensate each other within the integration domain allowed by 
the Gaussian WP (6.2). To clearly see these phases, it it convenient to change the integration 
variable according to p^ = + k^^ and then expand (with = —2k^k- + k^) 

\P\ = V-iQ + ky = VQ^ -2q-k-K^ 

^ ( q+k^ +q^k+ K"^ \ ^ /- . k+ k'] , ^ 

-Q(l + ^^-^-5^)-Q+V2,t+^-^. (6.5) 

This expansion requires a few words of explanation: among the subleading terms under the 
square root, the first one is of relative order (q ■ k)/Q'^ ~ l/{Qa) and hence it is much larger 
than the second one, which scales like K'^/Q^ ~ l/{Qa)'^. So, to be fully consistent, one should 
either push the expansion to the second order or neglect the last term, oc K^/Q^. However, our 
purpose here is merely to determine the dominant dependencies of the bulk excitation upon z, 
X- and x± at large x^. These are correctly encoded in the terms proportional to k-, k+ and 
respectively k^, as kept in the last approximate equality in Eq. (6.5). 

Specifically, using this approximation (6.5), we shall now successively perform the integra- 
tions over k-, k^ and k± in Eq. (6.1) (recall that we set p^ = q^ + k^). To that aim, we shall 
keep only the dominant, exponential, dependence upon k^ and replace k^ ^ (i.e. p^ — ?• q^) in 
the prefactors. The relevant integrals are then Gaussian and can be easily performed: 

dt- _/_^\ ^ _J_ 3,J_(f±:i^l , ,6.0) 
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Eq. (6.6) shows that the bulk excitation is itself a wave-packet which at time x"*" ~ is 
localized in the radial direction near z = with (recall that (7_|_ ~ 70") 



t 

7 



< 



(6.9) 



Eq. (6.7) shows that for a given z, the bulk WP is localized near x ■ = z/{2^/2'y) with an 
uncertainty C7_ ~ '7/7. Since moreover z ~ z^,, this implies 
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or 



t - X3 



2^ 



(6.10) 



This value of x„ is of course the maximal longitudinal extent consistent with Lorentz contraction. 
Finally, Eq. (6.8) together with Eq. (6.9) imply the following condition for the average position 
(x^) of the bulk WP in the transverse plane (note that z*/Q = t/q^): 



{xD ^a\ + {t/q^Y 



t 



(6.11) 



if t/gO>a2 



The second line shows that, for very large times, the bulk excitation expands in the transverse 
plane via diffusion. 

To summarize, the bulk excitation produced by the decaying system at time 3> a+ reads 
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Consider now the corresponding excitation induced by the probe. The respective boundary 
field is shown in Eq. (2.11), which immediately implies 



^'A(a;,z;T) 



(2^)4 

where (with C = C{2'Kfd+a-a\ 



(6.13) 



ipA{k', t) = C' exp 



iik_ 
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A+)^ aiik 



(6.14) 



As already mentioned, the central 4-momentum is taken to be space-like, A^ = A^ — 
2A+A_ > 0, and the widths of the WP are assumed to be large enough for the condition 
A:^ > to be obeyed by the typical modes A:'^ contributing to the integral in Eq. (6.13). Hence, 
the relevant expression for the bulk-to-boundary propagator is that given in the first line of 
Eq. (6.3). In practice, we shall take o"_ and a± to be so large that the respective momenta 
have only negligible spread: /c+ ~ A_|_ and k± ~ A_|_. This fixes the longitudinal and transverse 
resolution of the probe. As for the temporal resolution, this is controlled by the Gaussian in 
x+ in Eq. (2.11), which is centered at r with a width <C r. Accordingly, the central value 
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A_ of A;„ is not really needed and one can choose A„ = without loss of generality. Yet, the 
typical modes in the WP (6.14) will have a non-zero light-cone energy ~ 1/^+- So, a typical 
probe mode will have a virtuality = — 2A+A:_ with k- ~ l/<5'+- The condition that this 
virtuality ET^ remains positive (i.e. space-like) implies the constraint 

A+<a+Ai. (6.15) 

Eq. (6.15) plays the same role in the present context as (5.14) in the context of Sect. 5.2: it is an 
upper limit on the longitudinal resolution of the space-like probe. According to this equation, 
the best longitudinal resolution for a given A_l is achieved by choosing the largest possible value 
for At large time, it is convenient to let cj+ increase with r, like = eT, with e ^ 1. Then 
the constraint (6.15) becomes A+ <^ erA^. 

6.2 The three and four point functions 

The scalar versions of the three-point and four-point functions of interest, cf. Eqs. (4.9) and 
(5.10), read 

G(3)(g;r,A) = {dldA{T)d,), G(^)(g;r,A) ^ (Oj o^(^) o^(^) o^) , (6.16) 

where Oq and Oa are smeared versions of the 'dilaton' operator, as obtained after averaging 
over the respective (time-like or space-like) WP: 



Og = J d''x^g{x)Oix), Oa{t) = j d''xil^A{x]T)0{x). (6.17) 

As announced, the scalar toy model under consideration is characterized by cubic self- 
interactions. At tree-level, which is the relevant approximation in the context of SUGRA, 
this cubic interaction contributes to the three-point function in (6.16) via the Witten diagram 
shown in Fig. 6. a. We shall evaluate this diagram using the SUGRA Feynman rules along the 
Schwinger-Keldysh contour in the complex time plane (the 'closed time path'), as appropriate 
for computing quantum correlations in real time. Our use of the Schwinger-Keldysh formalism 
will however be quite minimal, so we shall not describe it here in any detail. (See e.g. [38, 39, 41] 
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for recent applications of this formalism to Witten diagrams in SUGRA.) For the diagram in 
Fig. 6. a, one finds 



dz d'^p d^/c 



As 



z^ (27r)4 (27r)4 



Dn{p, z) Duip + k, z) Dn{k, z) (t>q{p) </.*(p + k) ^^{h r), (6.18) 



with A3 denoting the strength of the cubic vertex and where we recall that y/—g = / z^ . Here 
the lower indices on the propagator refer to the branch ('1' or '2') of the Schwinger-Keldysh 
contour that the corresponding scalar field has been placed on. Dii{p,z) is the time-ordered 
bulk-to-boundary propagator, as already shown in Eq. (6.3). -Di2(p + k,z) is the Wightman 
(or 'cut') bulk-to-boundary propagator, which is non-zero only for TL momenta, in which 
case it can be obtained by taking the — )■ limit of the respective bulk-to-bulk propagator 
shown in Eq. (6.32) below. This yields an expression like the complex conjugate of the second 



line in Eq. (6.3), but with the replacement HV 
r(2) 



(2) 



-2U2 



-i(H2^'* +H2^^). For the present 



purposes, it is only which needs to be retained. (Indeed, we are interested in the behaviour 
of the integrand for relatively large values of z and the other function H2^^ would lead to strong 
oscillations in that regime; see below). To summarize, within the integrand of Eq. (6.18) we can 
replace 



Dii{p,z)Di2{p + k,z)Dii( 



i Z^TL {p, z) I^TL {p + k,z) DsL {k, z) 



(6.19) 



with DgL and Dtl denoting the SL and respectively TL components of the Feynman propa- 
gator in Eq. (6.3). The boundary WPs are shown in Eqs. (6.2) and (6.14). As discussed after 
Eq. (6.14) the probe (SL) WP is taken to be sharply peaked around and A_l, so that the re- 
spective integrations over and k± can be trivially performed. Regarding the TL propagators, 
expanding for large arguments and taking k '^p and ^ Q^, we have 



D-j:i,(jp,z) D^^{p + k,z) 



TTZ'Q' 



irQz 



exp 



\/2-fk- + 



k^ 



Then the p-integration is simply related to the normalization of the TL WP since 



d'^p 
(2^ 



d^p 
(2^ 



\Mp)\' 



d''x\Mx)\\ 



which is taken to be equal to unity. Putting everything together and defining = k^ 
and K'^ = A]_- 2A+(£_ + A_) ~ A]_ (recaU the constraint (6.15) on A+) , we obtain 



(6.20) 

(6.21) 

- A_ 



G(3)(g;A,r) = -iA3C<7+ I 



dz d^„ 7r{Qzf{Kz)\ 



V2^ 8 2 
exp \il-{T - V2-fz) 



die 



IZ 



( Ah 



V2y27 



+ \/27A 



(6.22) 



Now the ^--integration can be easily performed and we obtain 



G(3)(g;A,r) 



dz TT{Qzf (A^z) 



K2(Axz)e 



A+ 

472 



— 



25-i 



(6.23) 
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where it should be clear in the above that is not an independent variable, but simply deter- 
mined by Xj^ = \/2'yz. Also, we temporarily keep a non-zero value for A_, to explicitly show 
that this needs to be small, A_ < l/r. Note that the Gaussian restricts xj^ (the interaction time 
in the bulk) to values which are relatively close to r (the measurement time on the boundary): 
|x+ — t| < 0'+. Recalling our interpretation of x+ as the physical emission time (the time when 
the signal measured by the three-point function is actually emitted by the target, cf. Sect. 3.2), 
this confirms the fact that a space-like probe is a good measurement of the state of the system 
at the time of measurement. 

Before performing the final integration over z, let us open a small parenthesis and notice 
that Eq. (6.23) does not carry any information on the widths of the TL WP. This remark is most 
easily understood in a coordinate space calculation; since we are not interested in discriminating 
the internal structure of the TL WP, the probe resolutions have been assumed to be low on the 
scales set by the various widths in Eq. (6.12), that is, A+ <^ l/c-, A_l <C and » (J+. 

In view of this, we can replace (the modulus squared of) the bulk excitation (6.12) by its formal 
limit as obtained when all the widths approach to zero. This is a 4-dimensional (^-function 
(recah that \N\'^ = 1/ {ir'^a+a^a'])) : 

\%{x,z)\' = - 5(2)(^^)5(x+ - V2jz). (6.24) 

As discussed in Sect. 3.1, this simplification has already been used in the calculation of the 
backreaction. Eq. (6.24) represents a particle falling into AdSs in a highly boosted frame; 
indeed, the trajectory of such a particle is 

1 — V x+ x+ 

X- = — — z = —— — -x+~^=^ and a;_L = 0_l, 6.25 
l + v 4Y 7(1 + V) y/2j 

where the approximate equalities hold in the infinite momentum frame 7^1. 

The presence of the (5-functions in Eq. (6.24) greatly simplifies the calculation of the three- 
point function when using the coordinate space expression of (6.18). The probe excitation 
appears in the form ^^(2;+, a;+/47^, 0_l, z; r) with x-^ = ^/2^z. Then by using the coordinate- 
space version of the probe WP, Eq. (6.13), one easily recovers Eq. (6.23). 

Let us now close our parenthesis and return to do the z-integration in Eq. (6.23). We would 
like to argue that the prefactor in the integrand, to be succinctly denoted as /, is slowly varying 
within the range of the integration over z, that is, when x+ = \/2jz is changing from t to T + a+. 
For definiteness, we shall do this in the two limiting cases, A±z ^ 1 and A±z <^ 1. Recalling 
that z = x+/\/2j, we see that in the first case, namely A±z = A±x^/\/2j ^ 1, the prefactor 
behaves like / ~ x° with a = —2. (The actual value of a is not important for the argument, so 
long as this is not too large.) Then one has 

Sf a+f aa+ aa+ 

— ~ — — = ~ <C 1, (6.26) 

f f X+ T 

since we recall that cj+ ^ r. In the other case A±z S> 1, the dominant dependence is coming 
from the exponential falloff of the modified Bessel function and one has 

^ ~ ^ « 1, (6.27) 

/ 7 
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so long as A_|_ is not getting too large compared to 7/t. 

From now on, we shall work under the assumption A_|_ < 0(7/t), that we shall a posteriori 
check to be satisfied for all the situations of interest. Then the prefactor is slowly varying in the 
two limiting cases, as anticipated, and it keeps this property at all intermediate values. Hence, 
we can easily do the Gaussian integration to arrive at 

G(3)(,;A,r) = -a3V^^^^i^i^f^K2(A^..), (6.28) 

72;° 8 1 

where we have defined = Notice that in order to reach the above result we have 

neglected a factor originating from the phase of the space-like WP, namely 



exp 



A+ \ al /A+ ^ 



(6.29) 



This is correct so long as we impose that phase to be small, that is, A4, <^ ^'^ 1^ and A_ <C l/r. 
The first condition is automatically satisfied, as it is weaker than the constraint (6.15). Indeed, 
using cr+ = er with e ^ 1 and A_l < 7/r within Eq. (6.15), we deduce 

A+ < e-i'^jr. (6.30) 

The second condition A_ <C l/r is necessary to avoid strong oscillations and motivates us to 
choose A_ = 0. 

The presence of the modified Bessel function K2(A_|_z*) in Eq. (6.28), which originates from 
the boundary-to-bulk propagator for the space-like probe, effectively restricts the transverse 
momenta to A_|_ < 7/r, as anticipated. The result (6.28) is formally independent of A+, but 
clearly this is valid only for longitudinal momenta A_|_ obeying the constraint (6.30), which is 
the condition that the probe be space-like. The physical consequences of these constraints will 
be discussed in the next section. Also, notice that, in so far as the dominant behaviour upon 
A_|_z* is concerned, this result, Eq. (6.28), of the scalar toy theory is in fact consistent with the 
respective result of the backreaction, for the same type of bulk excitation (a falling particle) and 
the same, space-like, kinematics (compare to Eq. (3.23)). 

The Witten diagram contributing to the four-point function in the toy-model, scalar theory 
under consideration is shown in Fig. 6.b. In analogy to the three-point function, it can be 
estimated as 

G(^)(g;A,T) = A2 ^d^x^ I d\'^^,{x,z)^A{x,z;T)G{x-x',z,z')^l{x',z')^l{x',z';T) 

2 fdzdz' d'^p dV d'^k d'^k' -sirti), , , , ,,n 
= An / ^ — r 7 — TT . \, -, — TT -, — TT i^T^) ' ip + k- p - k) 
^ J z^ (27r)4 (27r)4 (27r)4 (27r)4 ^ ^ ^-^ f 1 

Gi2(p + k, z, z') Z?tl(p, z) D*^^{p', z') Dsj^ik, z) D*s^{k' , z') 

</'g(p)0$(p)V'A(fc;r)VA(^"';T) (6.31) 

where the only new ingredient is the Wightman bulk-to-bulk propagator: 

Gi2(£, z, z') = 7re(-£2) z'^z'^h{\t\z) h{W). (6.32) 
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The various bulk-to-boundary propagators visible in the integrand of Eq. (6.31) have entered 
the calculation as time-ordered (-Dii) or anti-time-ordered {D22) propagators in real time. 

Since A;, k' <C we shall shortly see that the dominant exponential dependence on p and 
p' in the propagators product cancels. Regarding the remaining dependence upon p and p' in 
the prefactors (as arising from the large argument expansion of the propagators), these are weak 
and hence we can simply replace p and p' by their central value q. Then one can integrate over 
p and p' to recover the normalization conditions for the TL wave-packets: 



d^p 
(2^ 



(6.33) 



(6.34) 



I WfWV' ^^""^'^^'^ ip + k-p'- k') Up) <t>W) ' 

At this stage, we managed to bring Eq. (6.31) into a factorized form: 

G(^)((/;A,t) =GL(<z;A,r)GR(g;A,T) = \GUq;^,r)\^ 
with the "left" part given by 

GL{q;A,T) = XsV^ I ^-^J^{\q + k\z)DTUQ,z)DsL{k,z)^A{k;r). (6.35) 

Prom this point on, the calculation is very similar to the one of the three-point function. For 
large arguments we can approximate 



32i\q + k\z)R^i\Qz 



1 2 



2ttQz 



exp 



-iz 



V2-fk^ + 



2^27/ 



+ 



(6.36) 



where we have used Eq. (6.5) and with the dots standing for a term proportional to exp(2iQz) 
which is neglected since it is varying rapidly. Taking again the SL WP to be sharply peaked we 
find (cf. the similarity with Eq. (6.22)) 



Gh[q; A, r) = z Vvr A3 C cr+ / ^ ^= — — K2{Kz) 



2tt 4 



exp 



iz 



( 



V 2^/27 



+ V27A, 



Performing the integrations over ^_ and z as usual we get 
G'L(<?;A,r) 



iX^Tl ^ — K2(A_L2:* 



. (6.37) 



(6.38) 



One should also include in the above a factor equal to that in Eq. (6.29), that is, the product of 
a phase factor times a Gaussian. However, when we construct the modulus squared according 
to Eq. (6.34), the respective phase factors mutually compensated^, while the Gaussian factors 
can be safely set to unity. We thus finally arrive at 



G<.)(,; A. .) = Al ,^ (%)! <^ [K,(A,..)1^ 



16 



(6.39) 



At this level, there is a small difference compared to the case of the three-point function: the two phases 
g±iA n (■j-ggg^jj jTjq (5,14)) automatically cancel out between the two insertions, $a and of the space-like WP, 
so the condition A_ < 1/r is not necessary anymore. 
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There is some model-dependence inherent in this formula, but this is harmless for our present 
purposes: the only information that we actually need is the dependence of the four-point function 
upon the dimensionless variable A_|_z* = A_|_r/\/27. This dependence is robustly predicted by 
Eq. (6.39) and could have been anticipated without explicitly performing the calculation, as 
we explain now. The four-point function defined in Eq. (6.16) represents the imaginary part 
of a forward scattering amplitude, which at the level of the SUGRA calculation is obtained by 
taking the 'cut' of the 4-leg Witten diagram depicted in Fig. 6.b. (This cut is manifest in our 
use of the Wightman prescription for the bulk-to-bulk propagator in Eq. (6.32).) In turn, this 
cut diagram is proportional to the modulus squared of the 3-leg diagram shown in Fig. 6. a and 
evaluated in Eq. (6.28). (This is merely the optical theorem adapted to the SUGRA context at 
hand.) We conclude that, at strong coupling, the four-point function G^'^\q; A,t) must depend 
upon A_lz* in the same way as the square \G^^\q; A,t)\'^ of the three-point function. This 
conclusion is indeed consistent with our previous results for G^'^\q;A,T) (see Eq. (6.39)) and, 
respectively, for G^^\q; A,t) (cf. Eq. (3.23) or (6.28)). Yet, this formal similarity between the 
three-point and the four-point functions should not be misleading: the physical content of these 
two quantities is very different, as it will be further discussed in the next subsection. 

6.3 Physical discussion 

The physical discussion to follow will only exploit those aspects of our above result, Eq. (6.39), 
for the four-point function at strong coupling which are firmly under control: its dependence 
upon the transverse momentum of the probe, which enters via the dimensionless variable A±z^, = 
A_|_r/\/27, and the upper limit (6.30) on the longitudinal momentum A4. of the probe. It is 
useful to recall that, in this boosted frame, the decaying system has an overall transverse size 
Ax± t/j and that the maximal longitudinal broadening permitted by causality and Lorentz 
contraction is Ax_ ~ t/j'^ (cf Fig. 1). We would like to check whether the system involves 
some substructures like partons with sizes much smaller than this maximal sizes. To that aim, 
one needs to estimate the four-point function for relatively high momenta A_l ^ 7/T and 
A- ^ 7^/^. (Such values are compatible with the space-likeness condition (6.30) provided A_l 
is chosen to be high enough.) But for such large values of the momenta, the four-point function 
(6.39) is exponentially suppressed, because K2(A_lz*) oc exp{— A_lz*} when A±z^ ^ 1. This 
simply means that there are no substructures in the decaying system at large times t ^ a^: the 
matter is uniformly distributed within a region with transverse size Ax± ~ r/7 and longitudinal 
size Ax- ~ T ll'^ ■ In particular, it exhibits maximal longitudinal broadening. 

It is amusing to notice that, at strong coupling, the four-point function and the three-point 
function are formally similar to each other (compare Eqs. (6.28) and (6.39)) — they both predict 
the exponential suppression of the respective correlation for transverse momenta Aj^ ^ ll^- 
However, as it should be clear from our previous analysis, this similarity is deceiving. The three- 
point function is independent of the coupling, so it looks quasi-homogeneous in the transverse 
plane (when probed with a low longitudinal resolution) even at weak coupling, where point -like 
partons are well known to exist. This is so since, by construction, the three-point function 
integrated over x_ represents the average energy per unit transverse area, which in this problem 
is homogeneous by symmetry. 

On the other hand, the four-point function has the potential to reveal small fluctuations in 
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the transverse plane, as manifest from the corresponding discussion at weak coupUng. So the 
corresponding exponential suppression for transverse momenta A_l ^ 7/t is a unambiguous 
proof of the absence of partons. 

It is furthermore interesting to re-express our results in terms of the typical momenta and 
virtualities of the quanta composing the decaying system at time r. By the uncertainty principle, 
a quantum distribution of matter which looks homogeneous (in a given event) over transverse 
distances 6x± ^ t/j involves Fourier modes with transverse momenta and virtualities k± ^ fj, ^ 
7/r and hence with longitudinal momenta fc+ ~ For comparison with the weak coupling 

discussion in Sect. 5, we note that the longitudinal momentum fraction x = k^/q+ carried by a 
typical mode at strong coupling is x ~ Xc{t) with 



This is independent of the coupling, unlike the corresponding weak-coupling result in Eq. (5.8). 
In fact, Eq. (6.40) looks more like the extrapolation of Eq. (5.8) to values of the coupling of 
order one, rather than to A — )• 00. This is consistent with the fact that the time dependence 
of X shown in Eq. (6.40) is the fastest one to be allowed by causality. Interestingly, the above 
estimate for Xc can also be written as Xc ~ fJ'/Q, which shows that, at strong coupling-*^^, the 
fraction of longitudinal momentum carried by the typical quanta in the decaying system is 
commensurable with the respective fraction of virtuality. This is in agreement with the picture 
of democratic parton branching, as put forward in Ref. [7] , in which the energy and the virtuality 
are quasi-democratically split among the daughter partons at any branching. In Appendix B, 
we shall present an alternative derivation of Eq. (6.40), which is in the spirit of the perturbative 
calculation for the fragmentation function in Sect. 5.1 — that is, it relies on the expression for 
the time-like anomalous dimension at strong coupling, as obtained in Ref. [23]. 

To summarize, at strong coupling and for sufficiently large times after the decay has been 
initiated, the decaying system occupies the maximal region in space and time which is allowed 
by causality and special relativity, that is, Aa;_|_ ~ T/7 and Ax_ ~ r/7^. In the center of mass 
frame of the virtual photon, the matter produced by its decay at time t ^ a is spread over the 
whole ball with radius r < t and its distribution within this ball is (quasi)homogeneous. The 
strongly coupled matter looks like a soft, smooth, jelly. 

Acknowledgments 

We would like to thank Simon Caron-Huot for an observation which triggered part of this work. 
One of us (E.I.) would like to thank the Physics Department at the Universidad Federal de Rio 
de Janeiro for hospitality and Prof. Nelson Braga for inciting discussions during the completion 
of this work. The work of E.I. is supported in part by the European Research Council under 
the Advanced Investigator Grant ERC-AD-267258. The work of A.H. M. is supported in part 
by the US Department of Energy. All figures were created with Jaxodraw [48] . 

^*For comparison, note that at weak coupling the virtuality fi is decreasing with time, via successive branchings, 
much faster than the longitudinal momentum fraction Xc, this can be checked using the general relation Eq. (5.7) 
together with Eq. (5.8) for Xc at weak coupling. 




(6.40) 



- 44 - 



A The Sb term for the faUing particle 

In Sect. 3.1 of [13] we have calculated the boundary 'energy density' £ generated by a particle 
falling in AdSs and with velocity v along the X3 axis. To this end, following [12] we have 
separated £ in two pieces £a and £3- We have erroneously stated that £b vanishes; this is true 
only when f = as we shall shortly see. Notice that, since in [13] we were in fact interested 
only in this particular case where v = 0, all subsequent calculations there were correct. 

We shall follow the notation of [13] with the only exception being the replacement Eq qq. 
In general, for an arbitrary 5D stress energy tensor tMN the term £b reads 

^""^l^J 0(i-i>'"(W) [(r - r)2(2too - 2*55 + - 3(x - x)\x - x)h,,\ , (A.l) 

where W = — (i — i)^ + — 't)^ + ■ For the falling particle under consideration we have 

7 'M 'N 

t'''' = qo {£j 6{x' - vt) 6^\x^) 5{z - ihf-^, (A.2) 

and by lowering indices with the metric in Eq. (3.2) we find that the non-zero components of 
interest are 

toQ = qQ(^^^ 5{df^ -vt)5^'^\x^)5{z-th), t33 = v\o, t55=too/7^. (A.3) 
Thus, the square bracket in the integrand of Eq. (A.l) becomes 

{r-rf{2too-2t55+tii)-3{x-xy{x-xytij = 3v^xl 6{x^ -vt) d'^^^x^) 6{z-t/j), (A.4) 

and using these 5- functions it is straightforward to integrate over x, x± and z to obtain 

POO 

d% / dtp 6{t^ - r2 - 2(t - vx3)t) . (A.5) 
Jo 



..2^2 roc 

£b 



2qo V 3 



vr 72 

As announced earlier, this contribution vanishes only when v = 0. Eq. (A.5) is the second term 
in Eq. (3.4) and the remaining part of the calculation leading to the final expression in Eq. (3.7) 
is given in the main body of the present paper. 

B Fragmentation function at strong coupling 

The time-like anomalous dimension at strong coupling reads [23] 



,0) = __^,_,„__^J. ,B.l) 

when j <C y/X and with jo = 2 — 2/\/A. The counterpart of Eq. (5.2) at strong coupling is 

^I?(x,QV/^')= [ ^Z)(j,l)e(^-2)"^(i/^)+^(^)''^(W) 
J 27ri 



x"^ ' 
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where the fragmentation function D{x^Q'^ / ^'^) is the special case of the four-point function 
computed in Sect. 6 for = fi^. Yet, it turns out that one cannot rely on Eq. (B.2) to 
recover the results of Sect. 6 in this special limit because of non-commutativity issues to be 
later explained. As in Sect. 5.1, we shall evaluate Eq. (B.2) using the saddle point method. The 
saddle point is located at 

For consistency this has to be much smaller than \/A, so that 

In^ = ln^-ln-<ln^, (B.4) 

/i fl X fi 

which requires x to be relatively close to n/Q; this condition is indeed satisfied, as we shall a 
posteriori check. Evaluating Eq. (B.2) at the saddle point, we get 

2nr 2, f ^ r o^2, Q'\ 2(QY I VA(lnxQ/^\ 
xI?(x,H~exp(^-^(j,-2) ln-^J=x ^-J exp 21nQ^ j' ^^-'^ 

This has a maximum at jg = 2 for which 

-c=f^V''^-^. (B.6) 



,Qy Q 

In terms of the formation time r = 2xcq+/ ■, we equivalently have (with tq = 2g+/(5^) 

xc^-, (B.7) 
r 

in agreement with Eq. (6.40). As A gets larger, the distribution Eq. (B.5) becomes more strongly 
peaked at x = Xc, but the limit A — )■ oo is subtle. Eq. (B.5) becomes the delta function 5{x — Xc), 
however, the limit A — )• oo and the j'-integral in Eq. (B.8) do not commute because the saddle 
point in Eq. (B.3) depends on A. Indeed, if we set A = oo first, the integrand of Eq. (B.8) 
becomes 

x'D{x, QV/^^) = / £ D{j, 1) , (B.8) 

for which there is no saddle point and the dependence of -D(j, 1) on j cannot be neglected. Still, 
one can see that the right-hand-side depends only on xQ/ ^ and decays rapidly as a function 
of this variable. In order to determine its functional form, one needs a direct evaluation of the 
fragmentation function in the context of SUGRA, as we did in Sect. 6.2 (see also Ref. [24]) . 
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